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•n ■ The consequences of assuming the third-generation Yukawa couplings are all large 

and comparable are studied in the context of the minimal supersymmetric extension of 
the standard model. General aspects of the RG evolution of the parameters, theoretical 
constraints needed to ensure proper electroweak symmetry breaking, and experimen- 
tal tal and cosmological bounds on low-energy parameters are presented. We also present 
(— I i complete and exact semi-analytic solutions to the 1-loop RG equations. Focusing on 
SU(5) or SO(IO) unification, we analyze the relationship between the top and bottom 
masses and the superspectrum, and the phenomenological implications of the GUT con- 
ditions on scalar masses. Future experimental measurements of the superspectrum and 
I of the strong coupling will distinguish between various GUT-scale scenarios. And if 
present experimental knowledge is to be accounted for most naturally, a particular set 
of predictions is singled out. 
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1 Introduction 



Unified theories are currently the most promising candidates for physics beyond the standard 
model. The marriage of force unification, namely grand unified theories (GUTs) or perhaps 
string theory, and spin unification, by which we mean a supersymmetry relating fermions 
and bosons, has been a fruitful and prolific area of research in the last few decades. Grand 
unification allows the understanding of the electroweak and strong forces as low-energy man- 
ifestations of a single microscopic force, in particular explaining the quantization and the 
assignments of electromagnetic charges for all known particles. The simplest GUTs [|l], 0, 
based on SU(5) or SO(IO) gauge groups, unify some [SU(5)] or all [SO(IO)] of the quarks 
and leptons in each generation. The unified matter multiplets neatly encompass the known 
standard-model matter particles — no new particles are needed, and no known particles are 
left out. The one exception is the right-handed neutrino, which must be included in SO(IO) 
unification. If neutrinos have masses, this potential embarrassment becomes a boon. In fact, 
SO(IO) not only favors typical ranges for their masses via the seesaw mechanism, but in 
specific models can also lead to detailed predictions about the flavor structure of the mass 
matrix. Moreover, SO(IO) beautifully incorporates both the Pati-Salam idea of leptons as 
the fourth color and an underlying symmetry between the left- and right-handed quarks and 
leptons 1^ . String theory aims to go beyond GUTs and unify all forces including gravity, per- 
haps without any adjustable parameters. Though still in their infancy, string models would 
presumably reproduce the successes of GUTs by either implying a grand unified theory as 
a "low-energy" consequence or by furnishing the appropriate boundary conditions to mimic 
grand unification predictions. GUTs, and even more so string theory, ambitiously span energy 
scales some 13 to 15 orders of magnitude above the highest scales at which the standard model 
has been directly tested. For compelling esthetic reasons such a span requires these theories to 
be supersymmetric: fermions and bosons, present in equal numbers, with mirror (and there- 
fore highly restricted) interactions. Supersymmetric theories are theoretically attractive on 
their own, but when wedded to unified theories they can yield quantitative phenomenological 
successes. Supersymmetry (SUSY) cannot be an exact symmetry of nature, but it must be 
approximately valid down to roughly the electroweak scale if such a low scale is to have a 
hope of being naturally embedded in a GUT- or string-scale theory The minimal can- 
didate model for unification is then given at low energies by the supersymmetrized standard 
model, which we will call the MSSMQ, having a squark for each quark, a slepton for each 
lepton, a gaugino for each gauge boson, and Higgsinos for the requisite two Higgs doublets 
and Hfi which give masses to the up- and down-type quarks and leptons (respectively). For 
a given effective theory below the unification scale, such as the MSSM, the renormalization 
group (RG) evolution determines the relationship between the physics of the unified theory 
and the physics observed at the electroweak scale. If one chooses as the effective theory the 
particular particle content of the MSSM, and embeds this minimal model in an SU(5) GUT, 

■^Here we use the term MSSM to refer only to the particle content and interactions of the minimal super- 
symmetric extension of the (essentially) experimentally extablished standard model, without any assumptions 
about the boundary conditions on its parameters at the GUT or string scales. 
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one arrives at a remarkably successful predicted relationship between the three low-energy 
gauge couplings ^. There are also various predictions for the quarks and leptons, some 
more robust than others. In particular, the bottom quark and tau lepton Yukawa couplings 
are in most cases predicted to be equal at the GUT scale 0; a successful prediction for the 
bottom mass at low energies, especially given a heavy top quark, can then be easily obtained 
with precisely the two Higgs doublets needed for the MSSM [0. The large Yukawa couphng 
of the heavy top quark can also trigger a correct breakdown of electroweak symmetry at low 
energies |jTl], |l^ . Thus the various pieces of the unification puzzle interlock tightly to produce 
a framework which we find both theoretically and phenomenologically compelling. 

Numerous authors have explored in detail the issues associated with gauge coupling uni- 
fication, in GUTs and lately also in string theories. The present paper takes a different path 
and seeks the consequences of Yukawa coupling unification 14 ] . We will only be concerned 
with the Yukawa couplings of the third generation, namely the top and bottom quarks and the 
tau lepton, for two reasons: first, because they are larger than the Yukawas of the lighter gen- 
erations, so it is natural to expect that they arise directly and from renormalizable operators 
and so are robustly predicted by the unified model, whereas the smaller couplings presumably 
arise from more complicated mechanisms and are thus more model-dependent; and second, 
because the third generation Yukawa couplings are the only ones big enough to appreciably 
influence the rest of the MSSM via the RG evolution. The focus of our research parallels these 



two motivations: we have seen in previous work |T5[ , as summarised below, that the top mass 



can be predicted from approximate or exact unification at the GUT scale; and we expand on 
our previous observations that the large bottom and tau Yukawa couplings which result 
from such unification qualitatively change the expected features of the MSSM at low energies. 

The assumption underlying this work is that, at the unification scale, either (I) A,- = A;, = 
At or at least (II) A^ = A?, ~ Aj, where \t,b,T are the Yukawa couplings to the appropriate 
Higgs doublet, and "~" means that the couplings differ by a factor of order one. When is 
this assumption valid? In the simplest SO(IO) scenario, in which both light Higgs doublets 
arise from a single SO(IO) multiplet, the tree-level Yukawa couplings are exactly equal at the 
GUT scale, as in assumption (I) |T^. Threshold corrections will typically lift this equality 
somewhat, and thereby can facilitate proper electroweak symmetry breaking, as we show 
below. In more involved SO(IO) models, the light Higgs doublets may come from mixtures 
of several SO(IO) multiplets. Nevertheless, we expect assumption (II) to often hold. In the 
simplest GUT scenarios based upon SU(5) the bottom and tau couplings are equal, but are 
unrelated to the top coupling. Most of the work on the MSSM has usually assumed that 
the top coupling was much larger than the other two, resulting in the observed hierarchy 
between the top and bottom quarks. From a GUT-scale model-building perspective it seems 
to us at least as natural a priori to assume that all three Yukawas are comparable, as in 
case (II); then the observed lightness of the bottom bottom and tau must result from the 
small vacuum expectation value (VEV) of the Higgs doublet to which they couple. Finally, 
the jury is still out on the predictions string theory makes for the Yukawa couplings. If the 
effective field theory which describes physics below the string scale is a GUT, then one of the 
above scenarios may hold [|l^]. Otherwise, there are still reasons to believe that the Yukawa 
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couplings for the third generation are approximately equal, at least in some string-inspired 
models. 

In addition to the above more theoretical motivations, there is also a phenomenological as- 
trophysical advantage to large tan/3, at least in SO(IO) models |T^, Various astrophysical 
and cosmological data, such as the neutrino solar flux deflcit and the density fluctuations at 
large scale observed by COBE, can be explained if the left-handed neutrinos acquire a mass 
via the seesaw mechanism from right-handed neutrino Majorana masses in the 10^° — 10^^ GeV 
range. A third-generation right-handed neutrino with such an intermediate mass can signif- 
icantly affect bottom-tau uniflcation (through the RG evolution of Yukawa couplings). If 
SO(10)-type boundary conditions = and \f = \% are assumed, it is difficult to re- 
produce the experimental value of mb/mr for small tan/5 (< 10 — 20). Interestingly, though, 
for larger values of tan (3, the intermediate-scale right handed neutrino does not significantly 
alter the successful prediction of rrib/m^. 

Motivated by the phenomenological successes of the unified MSSM and by the wide-ranging 
contexts displaying approximate or exact Yukawa unification, we analyze in this work various 
implications of a Yukawa-unified MSSM. In Sec. |^ we review the prediction of the top mass as 
a sensitive function not only of the GUT-scale boundary conditions and the bottom and tau 
masses, but also (perhaps surprisingly) of the superpartner masses. Reversing the argument, 
we find bounds on the superspectrum as functions of the top mass. Sec. | treats the related 
predictions for the radiative bottom quark decay 6 — 57 (and comments on r ^ /i7). We 
point out the importance of seemingly subdominant diagrams which are also enhanced by large 
tan (3. Satisfying the recent experimental bounds on this process places certain constraints on 
the superspectrum if a delicate fine-tuning is to be avoided. In Sec. || we outline the basic 
implications of A;, ~ At for electroweak symmetry breaking. Not only must the symmetry 
be broken radiatively without losing the SU(3)c x U(l)om gauge symmetries, but also a large 
hierarchy must be generated in the Higgs VEVs to account for the top-bottom mass hierarchy. 
Sec. ^ deals with the various options for generating this hierarchy. Two symmetries (PQ and 
R) can make this hierarchy more natural — and lead to a favored superspectrum — but there 
is always a necessary fine-tuning of at least one part in ~ 50 (~ mt/mb) somewhere in the 
Lagrangian [^, [1^. We return in Sec. ^to the problem of properly breaking the electroweak 
symmetry in the presence of the PQ and R symmetries. The various conditions which must 
be satisfied at appropriate scales to guarantee the proper spontaneous symmetry breaking are 
discussed in some detail in subsection In addition to the various individual mass-squared 
parameters, we examine the two fiat directions in the scalar potential (and the scales at which 
they may destabilize the vacuum) and the trilinear scalar couplings from A and fi terms. 
These can be important, even for the third generation, when there is some hierarchy between 
mz and the SUSY scale. We then turn to the general behavior of the soft scalar masses as they 
evolve down from the GUT scale, focusing in subsection |6.2| on the homogeneous part of the 
RG equations which dominates when the PQ and R symmetries are approximately valid. The 
favored scenarios with and without these symmetries are briefiy summarized in subsection |6.3| , 
without any assumptions about the soft scalar masses at the GUT scale. These assumptions 
are introduced in Sec. |^. The GUT-scale constraints on scalar masses in minimal SO(IO) 
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theories, and in SU(5) or nonminimal SO(IO) models, are presented in subsection using 



a common notation. The ramifications of these SO(IO)- or SU(5)-type GUT relations are 



explored in detail in subsections and ffTSl , respectively. These include the ease of obtaining 
proper electroweak symmetry breaking for different values of the parameters, the possible 
and probable superspectra, and the implications of lifting the PQ or R symmetries. We also 
study the effects of a right-handed neutrino mass below the GUT scale, briefly examine the 
consequences of universal scalar masses, and reconsider SU(5)-type boundary conditions as a 
perturbation on the SO(10)-type conditions. We then turn in Sec. ^]3|to the astrophysical and 
cosmological constraints on the Yukawa-unifled MSSM. We address both the electric neutrality 
of the lightest superpartner and its relic abundance. To estimate this abundance, we adapt 
previous analyses to the large tan/3 scenario, and in the case of a bino-like LSP classify 
the dominant operators contributing to LSP annihilation in order to clarify its suppression. 
We present our conclusions in Sec. 9. In particular, we summarize the phenomenological 
expectations from Yukawa uniflcation, comment on the most natural and therefore favored 
scenarios, and outline some directions for future investigation. In Appendix A we present 
the exact and complete semi-analytic solutions to the 1-loop RG equations for the MSSM 
with large Xt^b,T- They are semi-analytic in that they are given in terms of integrals over 
the dimensionless (gauge and Yukawa) couplings, which must be evaluated numerically or 
approximated analytically, as we show for several examples. Appendix B is devoted to a 
study of one of the potentially flat directions in the scalar potential of the MSSM, and to the 
scales at which it can impose a constraint on the scalar mass parameters. Finally, Appendix C 
justifles the approximation we have made in using the RG-improved tree-level scalar potential. 



2 Top mass: Prediction and Constraint 

At tree level, the observed masses of the third family fermions are related to their Yukawa 
couplings and to the VEVs of the up- and down-type Higgs doublets via: 

rut = XtVu = XfV sin (3 
rrij, = XbVo = XbV cos P 

rUr = XrVD = XrV cos (3 (1) 

where v = 174GeV and tan/3 = vu/vd = {Hfj)/{H^). The Yukawa couplings are in turn 
determined through the renormalization group evolution by the Yukawa couplings Xf^^ ^ at the 
GUT scale Mgut- And flnally, in the grand unifled theory these couplings are related to each 
other according to: 

A^ = A,« = Af = Ag, (2) 

or 

Xr = X'^ Xf = Xg . (3) 

In the minimal SO(IO) scenario A^/A^ = 1, while Higgs mixing or an SU(5) model could 
suggest that Af/Af is of order 1. The RG evolution requires as additional inputs the scale 
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of unification Mq and the unified value gc of the gauge couphngs, both of which are already 
fixed by gauge unification (but see Ref. for the treatment of a^). Thus the four low- 
energy observables mt, nib, nir and tan/3 and the two GUT-scale parameters Xg and Xf /Xf 
are related by three (RG) equations; fixing rrib and rrir from experiment leaves three equations 
in the four remaining variables mt, tan/5, Xg and Af/A^, yielding a single prediction for the 
top mass as a function of the angle (3. (In principle, of course, the Higgs doublets VEVs vu and 
vd and hence also (3 are predicted in terms of the GUT-scale parameters of the Higgs sector, 
but at this stage those parameters are completely unknown; we will return to them below.) 
If, as in minimal SO(IO), the initial ratio Xf/Xf is fixed, then the top mass and tan/? are 
predicted individually. We have previously presented a detailed study of this prediction fTS 



when Xf and A^ are split at most by threshold corrections; we will return to the more general 
case xf ~ A^ below. A crucial finding of that work (valid more generally for any Xf ~ A^) 
was related to the fact that, since the mass of the bottom quark results from a large coupling 
to a small VEV rather than a small coupling to a large VEV, any chiral symmetries protecting 
this mass should act on the VEV and not on the Yukawa coupling as in the usual case. Such 
approximate symmetries, discussed in detail below, are not respected by the parameters of the 
generic MSSM, which therefore exhibits large corrections to mh from 1-loop couplings to the 
other Higgs VEV (]15[; see also [^). In particular, we found that the usual suppression factor 
of l/lGvr^ in the leading 1-loop corrections is a priori completely undone by the enhancement 
~ tan (3 from the larger VEV of the up-type Higgs doublet. The two dominant contributions 
are given by the diagrams shown in Fig. 1. Keeping only these corrections (and therefore 
dropping similar but smaller corrections to m^-), the third generation mass relations now read 



mt 
rrib 
m^ 



Xt 
Xt 



XG,Xf/Xf V sin 13 



XcXf/Xf vcosp ( 1 + 



6mb 
mb 



(4) 



A^ XG,Xf/Xf vcosP 



where we have explicitly shown the dependence of the three Yukawa couplings at low energies 
on the two GUT-scale parameters. The exact form of the corrections to the bottom mass 
was given in our previous work. A useful approximation is given by 6mb/mb = (tan/5/50) 5b, 
where 



50 



/(x) = (1 — a; + xlna;)/(l — x) 
at the electroweak scale, aw = 



eff 



m? 



3 aG 



glf 



mt 



m, 



2X1 f 



efT, 



m 



eff. 



(5) 



and are the gluino mass and the term evaluated 
gl/'^T^ and as = gll^T^ are the SU(2) and SU(3)c coupling 
strengths respectively and rrigg = |(m? -|- m~) is the average of the squared masses of the 

SU(2)-singlet bottom squark and the SU(2)-doublet third generation squarks. [In the second 
(subdominant) term we used At ~ 2Mi/2 = 2f^?7i^ and approximated m? + m 



Q 



m? + m'~. 



Also, the expressions are considerably modified when one of the stop or sbottom eigenvalues 
becomes very small, but this will not be relevant for the cases we study.] We see that even if 
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Fig. 1. The leading (finite) 1-loop MSSM contributions to the b quark mass. 



Xf /Xf is fixed, for example in minimal SO(IO), we now have the additional unknown 6h, which 
precludes a separate determination of rrit and tan f3. In other words, if we don't know enough 
about the superspectrum to pin down 6b, we cannot fix Ac by comparing the prediction of 
with experiment, and hence we lose the independent prediction of the top mass. 

It will prove very helpful to understand the magnitude of the Sruh/mf, corrections, and in 
fact of all large tan f3 enhancements, from the point of view of symmetries. To this end we 
recall the Peccei-Quinn (PQ) and R symmetries introduced in Ref. |T5[. The PQ symmetry 
amounts to setting fi = , while the R symmetry requires the vanishing of gaugino masses 
Ml/2, of the SUSY-breaking trilinear scalar couplings Ai, and of the bilinear SUSY-breaking 
Higgs coupling B. If either symmetry were exact, then when the up- type Higgs acquired a 
VEV the down-type VEV would remain zero, so tan (3 would be infinite. (Of course we have 
in mind the usual scenario in which the Higgs mass matrix has a negative eigenvalue in the 
Hu direction only.) Also, down- type quarks and leptons such as the bottom and tau would 
be exactly massless to all orders. We will see that these symmetries are the key to making 
large tan (3 as technically natural as possible: just as the bottom mass can be made as light as 
needed by imposing the usual chiral symmetry, so the PQ and R symmetries can be imposed 
on the Lagrangian to varying degrees. Unfortunately, current LEP bounds set strong bounds 
on how natural the large tan/5 scenario can be |]T6|, as we discuss below. But PQ and R 
are still the key to alleviating as much as possible the need for fine-tuning, and are also useful 
for classifying the various superspectra and discussing their phenomenological consequences. 

If the symmetries are only approximately valid, one must specify at what scale this approx- 
imation holds: we will see that in very fine-tuned scenarios, the squark masses mo evaluated 
at the electroweak scale are much smaller than their values Mq at the GUT scale, so an ap- 
proximately symmetric GUT Lagrangian having /i ~ M1/2 <^ Mq could yield a spectrum at 
observable energies /z ~ M1/2 ~ having no observable symmetries. Thus for those cases 
we will distinguish between having PQ and R symmetries at all scales, and having them only 
at high scales. 



In Fig. 2, we present the results of a detailed 2-loop analysis, following Ref. 15|, of the 
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Fig. 2. The predictions of 2-loop RG evolution (plus 1-loop threshold corrections) for 
the running bottom quark mass (in GeV), the pole top quark mass (in GeV), and tan /?, 
as functions of the GUT-scale Yukawa coupling Xq. To compute threshold corrections, 
we assumed the preferred superspectrum of Eq. (p2|). The solid lines correspond to exact 
Yukawa unification, while the dashed lines indicate Xq = Xf = 2Xf as an example of 
approximate Yukawa unification. The values of rrib, after the finite corrections Srub are 
added, should fall within the shaded bands in order to agree with the bottom mass as 



extracted from QCD sum rules |15|. 



running MS bottom quark mass mf,{mf,), the top quark pole mass, and the approximate ratio 
tan j3 as functions of Xq within the perturbative regime. We chose two sample values of the 
QCD coupling as{mz), and considered either exact or approximate Yukawa unification, where 
the former is defined to be = while the latter is exemplified by A^/A^ = 2. (There is in 
addition a slight logarithmic dependence on the superspectrum, as calculated in Ref. |T^; for 
definiteness we have assumed the spectrum singled out below, namely all squarks and sleptons 
and the pseudoscalar Higgs at ~ 600 GeV, while the SU(2) gaugino and the Higgsinos are at 
~ rriz-) Also shown are the corresponding allowed ranges for m;,(m;,) extracted from the 
data on e~*"e~ — > bb using QCD sum rules ||23|. We use the values obtained in the analysis of 



Ref. [|T5|, which update those in Ref. [Q. The uncertainty on mf,(mb) is essentially theoretical, 
being dominated by our ignorance of 0{a^) corrections to QCD sum rules. In the absence of 
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Table 1: The consequences of Yukawa unfication using 2-loop RG evolution and 1-loop loga- 
rithmic threshold corrections from the preferred spectrum of Eq. (^2]) . For every value of the 
top quark pole mass we list the unified Yukawa coupling at the GUT scale, the minimim and 
maximum values of 6mb/mb needed to bring the bottom quark prediction into agreement with 
experimental data, and the corresponding values of the tan /3-independent quantity 6b. 

6b, we could read off an allowed range of Ag by requiring agreement between the theoretical 
and experimental values of mb] then both mt and tan (3 could be predicted within some range. 
Instead, we can only determine the top mass (and tan/3) as functions of 6b- 

Turning the argument around, for a given top mass we can calculate the amount of cor- 
rection needed to bring the bottom mass prediction into agreement with experiment. 
We can then remove the tan /3 dependence, leaving only the spectrum-dependent quantity 6b- 
Table 1 displays these minimal and maximal allowed values of and 6b- (Actually, these 

bounds on and 6b themselves depend on the spectrum due to threshold corrections, but 
this dependence is a weak logarithmic one; typically, the logarithmic variation in mt is at 
most a few GeV for the more interesting higher values of mt- To obtain precise predictions, 
though, all thresholds corrections should be included using a definite superspectrum.) We 
learn that a positive 6b must be quite small, while a negative 6b may be large enough in mag- 
nitude to bring high predictions of m^ back into agreement with experiment. For example, 
when as{mz) = .115 and Xf = Xf, superpartner spectra for which \6b\ ^ 5% allow a light 
top quark, whereas spectra for which \6b\ ^ 16% mandate a light top, where by light we mean 
mt < 160 GeV. Conversely, when Xf = 2Xf, a superspectrum for which \6b\ ^ 15% mandates 
a heavy top and favors a small as{mz)- Fig. 3 translates this information into constraints 
on mfr, (the mass of the gaugino superpartner of the W) and fi at the electroweak scale. 
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Fig. 3. The values of ^ and the wino mass, normahzed to a typical squark mass (specif- 
ically uiq = ^(m? + w-^))) which allow proper bottom-tau unification for large tan/3. 
The shaded and hatched regions correspond to exact and approximate Yukawa-unified 
scenarios, respectively, while the solid and dashed lines represent the lower bounds im- 
posed by LEP chargino limits assuming ttiq is 600 (the preferred value) or 200 GeV, 
respectively. 



normalized to a typical squark mass mo (taken specifically to be rries), for various values of 
the top quark mass, agimz) and Af/A^. Also shown are the lower bounds on m^/mo and 
n/rrio imposed by LEP chargino limits, for mo = 600 GeV (solid lines) or mo = 200 GeV 
(dashed lines). Evidently the large, nonlogarithmic threshold correction 6b is actually of great 
use: since, unlike the typical logarithmic corrections, it is very sensitive to the superspectrum, 
we can use experimental measurements of the top mass to learn about the hierarchies in the 
superspectrum. And from a model-building viewpoint, we can exclude those theories in which 
the gaugino mass and /i parameter do not fall into the allowed ranges shown in Fig. 3. 

^ we 



To reemphasize the importance of the corrections to the bottom mass (see also p5 
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have studied the consequences of bottom-tau unification for arbitrary Af/A^, imposing only 
that Ac = Af < 2 to ensure the vahdity of a perturbative expansion. (A larger Af would 
result in a weak-scale top mass very close to the fixed-point value, regardless of bottom-tau 
unification. This fixed-point value depends on A^ and A^, as well as on a^.) In Fig. 4 we show 
the predictions for tan/? and 'm^°^'^, fixing the values of ^^(m^) and mh{mh) and integrating 
over all possible superpartner and pseudoscalar Higgs masses between niz and 1 TeV. The 
black regions correspond to 5b = 0, while in the gray regions —25% < 6b < 20%. The effects 
of finite 6b are striking. 



3 6 ^ S7 

There is another immediate phenomenological implication of the large 1-loop corrections in the 



tan P ^ 1 framework |T^, ^ : the same diagrams, but with a photon attached in all possible 
ways and with a fiavor- changing vertex as shown in Fig. 5, contribute to the bottom quark 
radiative decay b sj. These contributions, which for small tan/5 are typically somewhat 
smaller than or comparable to the 2-Higgs standard model contribution, are parametrically 
enhanced by a factor of tan (3 ~ 50 in the amplitude. But the CLEO bound ||28| on the 
inclusive branching ratio BR(6 — > 57) < 4.2 x 10~^ (at 95% CL) is already roughly saturated 
by the 2-Higgs standard model amplitude even if the charged Higgs is rather heavy (and is in 
fact oversaturated with a light charged Higgs). Hence the large tan/? contribution must not 
be too large. 



For two reasons p9|, we will focus our attention on the Yukawa-coupled ("primary") 
chargino-exchange diagram of Fig. 5(a) rather than the gluino-exchange diagram of Fig. 5(b) 
or the gauge-coupled ("secondary") chargino-exchange diagram of Fig. 5(c). First, throughout 
the relevant regions of parameter space, the primary chargino-exchange amplitude exceeds 
or approximately equals the other amplitudes in magnitude. Second, it is predominantly 
determined by the third generation, namely by stop exchange, and hence its magnitude is 
fixed by the Kobayashi-Maskawa quark-mixing matrix element Vts — Vcb', in contrast, the 
other two diagrams arise from squark mixing between the second and third generation, and 
therefore depend on an independent mixing angle, which we shall call V23. The primary 
chargino diagram also depends on the A parameter which mixes the SU(2)-doublet and -singlet 
stop squarks, but as we noted above, the RG equations typically fix A at low energies almost 
entirely in terms of the gaugino mass (independent of the GUT-scale value of A). So the 
sign and magnitude of this diagram is completely calculable M, ITB| in terms of the same 



parameters which enter 6b', we have found that, when 6b < 0, as must be the case for any 
sizeable l^^l, the chargino exchange amplitude A~- interferes constructively with the standard 
model and charged-Higgs amplitudes (^sm and Ah-, respectively). Hence there can be no 
cancellations between these, and the constraint on A~- is more severe. On the other hand, 

A 

the new angle V23, which determines the gluino and secondary chargino amplitudes A~ and 
A'~_, gets not only a contribution ~ Vcb through the RG evolution, but also one from the 
fiavor structure at the GUT scale. Since we would like to remain as model-independent as 
possible, we will not make any assumptions about this structure, and thus V23 will not be 
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Fig. 4. The ranges of top quark pole mass and tan /3 allowed by bottom-tau unification 
at MguT; for different values of the strong coupling and of the running bottom quark 
mass. We scan over all possible superspectra between mz and 1 TeV, allow arbitrary 
Xf/Xf, and restrict Ac = < 2 to ensure perturbativity. The effects of the 1-loop 
corrections 5b are evident by comparing the black regions, for which Sb = 0, with the 
gray ones, in which —25% < 6b < 20%. These predictions are expected to be accurate 
to within a few GeV. The dashed horizontal line and hatched band are the top mass 
central value and an estimate of its uncertainty, respectively, as recently announced by 



the CDF collaboration |26|. 



determined. By tuning the flavor parameters and thereby V23 one could cancel the various 
large tan f3 contributions to 6 — S7 against each other and avoid any bounds from this process. 
However, one should also be careful about other phenomenological implications of this new 
source of flavor violation. Since we are dealing with a grand unified theory, there is also a 
leptonic counterpart of the new mixing angle V23, and this gives relevant contributions to the 
rate for r — nj. A more detailed analysis of the potential for cancellations and of r(r — > /X7) 



is presented elsewhere PP[. Our approach here will be to take into account only the sizeable 
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Fig. 5. The primary chargino, gluino, and secondary chargino diagrams which contribute 
to radiative bottom quark decay and are enhanced by large tan /3. Note that while the 
amplitude in (a) is completely predicted by third-generation parameters and Vet, the 
other two depend sensitively on intragenerational squark mixing parameters and may 
be equally important. 



and calculable primary chargino contribution: in those parameter ranges for which it is small, 
the other amplitudes are typically also small and there is no conflict with experiment, while 
in those ranges for which ^sm + Ah- + A.~- exceeds the experimental bounds, the other 
diagrams must be tuned to cancel against these amplitudes. Our goal is to study when and 
to what degree such a cancellation is needed. 

We will use the following expressions |^ to calculate the limits on the MSSM parameters 
for large tan/3: 



BR(6 ^ 57) ~ BR(6 ^ cez/) 



I{mc/mb) [1 - {2/3TT)asimb)fQCDirnc/mb)] 



(6) 



where BR(6 cez/) = 0.107, r] = asijnz) / Oisijrii,) = 0.520 [using as{mz) = 0.12 as a reference 
value], C = —0.191, I{mc/mb) = 0.45, fQcuij^clf^h) = 2.41. The dimensionless amplitude is 
given by = Asm + Ah- + A~- where 



Ah- 
Az- 



3m^.(i)/ 
2m2 



mr 



m 



w , 



m 



tan [3 



m~ 



i2 



/i 



/ 2 

2 



(7) 
(8) 

(9) 
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and the various functions are 



= 36(x-l)3 +6(^^^"" 
^ ^ 6(x-l)2 + 3(x-l)3''''' ^''^ 



6(1 -x; 



3 



5 - 12x + Tx^ + 2x(2 - 3x) In x . (12) 



In Eq. (I) we assume that the two stop mass eigenvalues are roughly degenerate. (This is 
in particular a good approximation for the interesting situation in which the stop is heavier 
then the top and the diagonal soft masses are almost degenerate.) Notice that the crucial 
off-diagonal ijiiL mixing has been factored out in A~-. We have kept the exact dependence 
on the chargino mass matrix: 

U*f mwV^^in(3\^_^ ^ ( m~- \ ^^^^ 



y2cos/3 /i ; I 



m 



X2 



We can now determine how large A~~ can be without requiring some cancellation to 
avoid conflict with experiment. We find that if A~- = then the 2-Higgs standard model 
exactly saturates the experimental bound for itlh- = 1 TeV and rrit = 170 GeV; for those 
values Ah- — .15^sm- Assuming some theoretical uncertainty allows the charged Higgs to 
be significantly lighter: a 30% uncertainty in the branching ratio would allow a charged Higgs 
of 350 GeV {Ah- — -S^sm), while with a 50% uncertainty the charged Higgs can be as low 
as 250 GeV {Ah- — .75^sm)- For the ranges of top quark mass we are considering, the rrit 
dependence is much smaller than such theoretical uncertainties. Now, if we add constructively 
a supersymmetric contribution equal to 50% of the standard-model amplitude, the branching 
ratio is 30% above the experimental bound without including any charged-Higgs amplitude, 
or 50% above the experimental bound if we include a 25% charged-Higgs contribution. (These 
values correspond to taking the soft SUSY-breaking scalar masses to be roughly degenerate at 
~ 700 GeV while the fi and wino mass are at ~ mz, which we argue is the favored scenario.) 
The charged Higgs mass cannot be made too large without fine-tuning the Z mass, so a 
25% charged-Higgs contribution is a reasonable lower bound. Thus, to be conservative, we 
conclude that there must be some cancellations whenever A~- > .5^sm, either from the other 

A 

diagrams [Figs. 5(b) and 5(c)] due to tuning of the flavor physics or from cancellations in 
the A parameter due to tuning of the GUT values of A and the gaugino mass. We define a 
measure of the need for cancellations asQ: 



eb_s-y = mm 



SM 



1 . (14) 



Since A~_ — -5 Asm requires no more tuning than A~_ < .5.4sm, we define et^s-y so it saturates when 
A~_ = .5.4sM- For similar reasons we will define es = min(i3/m^, 1) and es = m.m{m'^ / m\, 1). 
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For example, if the chargino amplitude is 10 times greater than the saturating amplitude 
.5^smj then some other contribution must be adjusted with a precision of 0.1 to cancel against 
it. In fact, we will never need the exact definition of eb^s-y', all that we will require of it is that 
it be small whenever the chargino amplitude is considerably too big. 



4 Electroweak Symmetry Breaking 

Thus far, the existence of the electroweak-breaking Higgs VEVs vu and vd has been simply 
assumed. We must study the generation of those VEVs in order to understand how compatible 
is the idea of Yukawa unification with radiative electroweak breaking, how the SO(IO) or SU(5) 
symmetry constraints on the soft SUSY-breaking mass parameters affect this breaking, how 
the top-bottom mass hierarchy is obtained, and how natural or unnatural is such a scenario. 
This study will reveal favored ranges for the Yukawa couplings and the soft SUSY-breaking 
masses, and highlight the central role of the PQ and R symmetries in discussing electroweak 
breaking for Yukawa-unified models. 

Electroweak symmetry breaking is governed by the scalar potential for the neutral com- 
ponents of the two Higgs fields Hjj and Hd, which at tree level is of the form 



\Hu\^ - iHnff (15) 



gl + gi (\rr 12 



where g2 and gi are respectively the SU (2) and hypercharge gauge couplings. The parameters 
''^u D — l^lj D 1^'^ contain the soft-breaking Higgs masses as well as the fi parameter from 
the super potential, while B is the soft-breaking bilinear mass parameter. The conditions for 
proper electroweak breaking are well-known: 

ml + ml > 2\i^B\ (16) 

ensures that the potential is bounded from below, and 

mlml < i/B"^ (17) 

guarantees the existence of a minimum away from the origin and so breaks the symmetry. 
In practice, since \^B\ will always be much less than or at most comparable to and 
|m|)|, we can reduce these requirements to m\ = mfj + mj^ > (using the expression for the 
pseudoscalar Higgs mass) and mfj < (noting that large tan P means that the up-type Higgs 
gets the large VEV). 

It is useful to interpret the two above conditions for proper symmetry breaking in terms 
of the usual custodial symmetry exchanging up- and down-type particles. In practice we 
need < < m|,, which represents a substantial violation of this symmetry. In the usual 
scenario, with the initial condition Xt ^ Xb^r, this custodial violation is provided by the Yukawa 
couplings themselves. The large top Yukawa coupling drives negative the soft-breaking mass 
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parameter of the "up-type" Higgs Hu to which it couples, while the other Higgs is largely 
unaffected. Hence, even with custodially symmetric soft masses at Mq, it is very natural to 
obtain 0, [l^, at the electroweak scale, < while the other Higgs mass is sufficiently 
positive to ensure mfj + m|) > 0. On the other hand, with unified or approximately- unified 
Yukawa couplings {Xf/Xf ~ 1), it is hard to see why the two Higgs mass parameters should 
run differently, so the natural explanation for mfj < < mfj + mj^ is lost. In fact, with the 
boundary condition Xf = Xf, the only sources of custodial breaking in the couplings are the 
hypercharges and the absence of a right-handed neutrino (but presence of a right-handed tau). 
These furnish only a tiny splitting even when integrating from the GUT to the electroweak 
scales. If such a splitting between mfj and mj^ is to be ~ m| then the soft-breaking masses 
themselves must be considerably bigger than m|. To break the symmetry more naturally, 
custodial breaking must be enhanced. To this end, either mfj can be split by various means 
from mjj) already at the GUT scale, or one may relax the requirement Xf = Xf, which in any 
case is bound to be modified — either by a little, due to GUT-scale threshold corrections, or 
by a lot, in the case of mixed-SO(lO) or SU(5) models. 

Furthermore, in the usual scenario, if all mass parameters in the scalar potential are 
comparable then so are the VEVs of the two Higgs doublets; but comparable Higgs VEVs 
are perfectly compatible in the usual scenario with the experimental hierarchy mj ^ mb since 
that is furnished by the assumed hierarchy of Yukawa couplings ^ A;,. By contrast, in 
a unified scenario it is the ratio of VEVs which must be large. Upon minimizing the scalar 
potential Vq, one easily obtains (for large tan/3) 



vn fj,B jiB 



tan 13 Vu m\ + mjj, m^ ' 



(18) 



as well as m\ = —2mfj (which sets the scale). A large VEV hierarchy requires a small coupling 
between the two Higgs doublets, namely jiB <^ mfj + m'jy = m\, so that an expectation value 
in one is only weakly fed into the other. But this small Higgs coupling, as we now discuss, 
implies [0, |l6l a necessary degree of cancellation between some parameters at the GUT scale. 



5 Generating a Hierarchy 

From Eq. (|TB]) it is clear that to generate the hierarchy of VEVs (tan (3 ~ 50) necessary for 
the top-bottom (and top-tau) mass ratio, we need a small fiB or a large m\. This is difficult: 
on one hand, m\ = mfj + m|, and m| = —2mfj run quite similarly and are tightly coupled 
through the RG evolution, so it is difficult to make m\ much larger than the electroweak scale; 
on the other hand, fi cannot be much below the Z mass since otherwise a light chargino (or 
neutralino) would have been detected at LEP, and the RG equations imply that B is naturally 
at least as large as m^ — which again cannot be much below the Z mass without producing a 
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light chargino or neutralino. To make these arguments concrete, we can rewrite Eq. ([18| ) as 

1 jjB l^^yY ^ "^1 



50 m\ m| m\ 

> 0.9 = 



(19) 



The most natural scenario would display approximate PQ or R symmetries, with either fi -C 
mz or ~ 5 ~ A ^ mz or both, while all other masses would be around mz- When 
tan/? ^ 1, such small /i or result in a chargino much lighter than niz- Alas, this scenario 
is now experimentally ruled out: specifically, LEP |3l|] places a lower bound of roughly | mz 



on the lightest chargino mass. A lower bound of | mz on the lightest eigenvalue of the chargino 
mass matrix in Eq. (13) translates into a bound on the Lagrangian parameters of 

- (|mz)1 L^- (imz)1 >2mlr{lmzf. (20) 



m 



w 



Subject to this bound, the prefactor ^my^/m\ in Eq. (|TU|) is minimized when /i 
yielding iJLm^/m\ > 1/ A + 2~^/'^mw /mz — .88. Therefore, the parameters at the GUT scale 
must be adjusted so that, at the electroweak scale, either B is much lighter than its natural 
minimum value ~ fn^yi or m| is much less than its natural value ~ m\. (Note that in the first 
case tan (3 is the quantity that is tuned to be small, while in the second case m| is tuned to be 
small.) We quantify these tunings by and e^, respectively [see the footnote for Eq. (0)]. 
Thus the bounds from LEP imply that some GUT parameters must be adjusted to cancel 
with an accuracy of at least 



eB^z ^ — — ^ ~ — ; saturated when/i ~ m^ ~ mz (21) 



1 

tan (3 50 

and the least amount of tuning is required when n and are roughly as light as they can 
be, namely both near the Z mass. 

Eq. ( |T^ ) by itself does not distinguish whether B or m| should be tuned small; that is 
decided by e^^^^, which requires considerable tuning of the fiavor sector or of the A parameter 
at the GUT scale if the superpartners and the pseudoscalar Higgs are near mz- Table 2 
sketches four possible scenarios, distinguished by whether PQ and R are good symmetries (at 
low energies — rrij are the sfermion masses at the electroweak scale). Note that es^z favor the 
PQ- and R-symmetric case or the asymmetric case having all masses near m^, but that e^^^^ 
favors the former over the latter. Thus the most natural scenario as measured by these three 
criteria is the maximally symmetric one: a Lagrangian (|]16[; see also the first reference in 



3^ ) which is maximally PQ- and R-symmetric while obeying the LEP bounds and keeping 

(22) 



the superpartners as light as possible, that is, 

J mA ~ rriQ ~ ^ taxi [3 mz (~ 600 GeV) 



^ A B m^ ~ ~ mz 



where by mo we mean the typical mass of the squarks and sleptons evaluated at the electroweak 
scale. (To reiterate: we chose ~ ~ ^mg because this is the most natural case allowed by 
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no R 



R 



mo 



1 



mz ~ y mo 
rrii mo 



/i 



ymo 
mz ~ ^ 



mi 



ymo 
mo 



1 



es ~ 1 



PQ 



either 



mo 
mo 



m^ ~ f mo 
m ~ mo 

or 

niz ~ mo 
rrii ~ mo 



/i 



mo 



mz ~ 7 mo 
mj ~ mo 



_i_ 

50 
1 

J_ 
50 



50 
J_ 
50 



50 
1 
1 
7 



no PQ 



Table 2: Typical scenarios, and consequent fine-tunings, for generating the Higgs VEV hier- 
archy vd ^ l/50vu with or without the PQ and R symmetries. 



LEP — if ijl or m^ were smaller, mz would require further tuning to make it sufficiently light, 
while if they were much larger we would lose the advantages of the PQ and R symmetries 

and the tuning would again be exacerbated.) This spectrum implies a small correction to the 
bottom mass, 6b ~ 5%, and hence (from Table 1) a heavy top and preferably a somewhat low 
value of as{mz)- 



6 Correct Symmetry Breaking 
6.1 The General Problem 

We return now to the question of how the electroweak symmetry may be correctly broken, 
while preserving the SU(3)c x U(l)em gauge symmetries. In principle, what needs to be done 
is to study the effective potential V^.^ for field values ms < (f) < Mq, where ms denotes 
collectively the soft SUSY-breaking masses and (j) is the set of scalar fields in the MSSM. 
In practice, since we are dealing with a pcrturbative theory, we need only consider the RG- 
improved tree level potential Vo(A) renormalized at scales A between nig and Mqut- First, 
we must make sure that the potential is stable for A ^ ms, i.e. that Vo(A) is bounded from 
below at high scales. If this were not the case, the scale of gauge symmetry breaking would 
be ~ A ^ nis, which is phenomenologically unacceptable, and possibly even SU(3)c x U(l)em 
would be broken. Second, we need to guarantee that at scales A ~ ms the minimum of Vo(A) 
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is such that Hu and Hd, and no other fields, acquire nonvanishing expectation values. This 
amounts to imposing some positivity constraints on the sfermion mass-squared parameters, 
which we discuss below. Finally, since we are interested in large rather than small tan (3, the 
instability should arise in the Hu direction, while the small Hd VEV is generated through the 
mixing mass parameter Bfi. 

An essentially technical comment is in order here. We will be mainly concerned with the 
parameters related to the third family and to the Higgs sector. This is because, in the limit 
in which fiavor mixings and the Yukawa couplings of the two light families are neglected, the 
SUSY-breaking masses of these families have a numerically small impact on the RG evolution 
of the parameters of the third generation and the Higgs sector. The only effect @] is via 
the hypercharge D-term S (see Appendix A) , which is small because of the small hypercharge 
gauge coupling. Moreover the effects of S are completely determined by its GUT scale value 
Sg, since it renormalizes multiplicatively. Indeed with SU(5) or SO(IO) boundary conditions 
on scalar masses, Sg itself is completely specified by the soft masses of and Hd (since these 
are the only light incomplete matter SU(5) multiplets). 

The relevant parameters then consist of the seven soft-breaking scalar masses fifj = m\j — 
/i^, /if) = m\) — fi^, mf, rn^, m~, m| and m~, the three trilinear soft-breaking parameters 

b T-, the single bilinear soft-breaking parameter B, the /i term in the superpotential, and the 
three gaugino masses. Their 1-loop RG equations are given for reference in Appendix A. 

Let us now discuss in more detail the constraints which these parameters need to satisfy. 
We begin with those which must be satisfied at scales A ^ ms- It is well known that the 
MSSM, like any generic SUSY model, possesses a host of "approximately fiat" directions in 
the space of scalar fields 0. By "approximately fiat" we mean that the potential, at the 
renormalizable level, is only quadratic in those directions. In general, though, irrelevant 
operators suppressed by inverse powers of a large mass M such as Mqut or Mpianck can 



give an additional stabilizing contribution [^. To be conservative, we will always assume 
the superpotential contains an appropriate operator of the form Then along any 

such direction parametrized by a field (pa, the potential is essentially (see below) given by 
m'^\(j)a\'^ + \4>af/M^y where is equal to a sum of squared masses. Regardless of the sign 
of rria, there is no minimum for (pa ^ y/rriaM = Ahigh (~ 10^ GeV for M = Mqut), so 
parameters normalized at scales A > Ahigh can never yield an unwanted minimum. For 
(pa < Ahigh, the potential in dominated by the quadratic piece, though there may be a scale 
Alow below which a linear term may again stabilize the potential. In the absence of a linear 
term, the lowest scale of interest is Alow ~ ^s, at which the superpartners are integrated 
out. If rria > were to become negative at a critical scale Alow <^ A^ -C Ahigh, dimensional 
transmutation would take place: the VEV of (pa would be fixed by the one loop correction 
to the effective potential Vi(A) to be of order Ac (times a coupling constant). (Notice that, in 
the absence of the irrelevant operator, if were to be negative already at the GUT scale, then 
we would clearly be expanding around the wrong vacuum in the GUT theory.) In order to get 
acceptable low-energy physics we have then to impose m^(A) > 0, for all approximately-flat 
directions a and for all scales A between the Ahigh and Alow relevant to that (pa- 

When we restrict our attention to the flelds of the third family and the Higgs sector, there 
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are only two such flat directions: 



1. {Hu) = {Hd) = (pi, with all other fields at zero; and 



2. {H^) = 02, {L) 
zero 



+ 02Ai/Afe)^/^ {Q) = {be) = (02/u/Afe)^/^ with all other fields at 



36 



The color and isospin orientations are determined by imposing vanishing D and F terms. 
Along 01 the potential is purely quadratic, Vo(0i) = mi|0ip with ml = + w?^ — 2|i?/i|. 
The stability constraint has already been given in Eq. (|16|) . This constraint should be satisfied 
between Ahigh ~ 10^ GeV (to be conservative) and Alow ~ "^5- Along direction 02 there is 
also a linear term: 



K)(02) = (/i^ + mf)|02p + (m| + m?5 + mf) 



2| J, |2 I 2 
m2|02| +"^3 



/i02 



A. 



A, 



(23) 



For this fiat direction, the dominant stabilizing term at high scales is the left-handed neutrino 
mass operator (HuL)'^ /Mn, where Mn is the right-handed neutrino mass. Indeed the effect 
of this operator can be important down to Ahigh ~ 10^ GeV, since could be as low as 
10^^ — 10^'^ GeV. At low scales, the linear term will stabilize the potential (provided m| > 0, 
which will always be the case). Therefore, as we show in Appendix B (see also Ref. [0), the 
02 fiat direction can only pose a problem at scales above Alow ~ (2vr/a)yu/Ab ~ 10'* — 10^ GeV. 
So we need to impose m2{A) > at least for all A between Ahigh ~ 10^ GeV and Alow ~ 
10^ GeV. 

A general scan of the values of mi2(A) at all intermediate scales would be numerically 
arduous. Fortunately, with only minor assumptions on the initial parameters, mKA) and 
m|(A) decrease essentially monotonically with A. Imposing positivity just at low A then 
guarantees the absence of unwanted vacua at all scales. Consider for instance the PQ- and R- 
symmetric limit of the RG equations in Appendix A. Monotonicity of ml ^ is clearly satisfied 
when Xt^b,T are positive throughout the running. In turn this condition is satisfied when the 
Xi start out positive and of comparable magnitudes (check for instance the entries in the 
matrix Ti in Appendix A, whose behavior in monotonic in A). In most interesting cases, 
the necessary positivity of masses at low energy will imply positive Xi at the GUT scale 
[for instance in minimal SO(IO)]. Introducing a finite /i does not alter the conclusions, as 
long as R symmetry is preserved. For small /i and large gaugino masses, the situation is 
also unchanged: in the first stage of the running, their contribution to m\ 2 is positive, but 
very small; however it soon becomes negative and its absolute value increases monotonically 
when A is lowered, so again checking positivity of m\ 2 at low scales suffices. Finally, when 
both PQ and R are broken, the above discussion applies straightforwardly to m\, but not 
to m\ due to the additional inhomogeneous piece B^. For this situation we have explicitly 
verified monotonicity for a wide range of initial parameters. We thus conclude that, quite 
generally, the imposition of the constraints at a low scale is sufficient to ensure their validity 
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throughout the RG evolution. Our analysis is thereby considerably simplified: we need only 
impose (A = ms ^ mz) — m\ > and 



to avoid an instability in the 0i and 02 directions. 

We note in passing that the constraints from fiat directions involving also the fields in 
the first two families are not a problem for us. This is because their always involve the 
soft masses for these fields, which are for us essentially arbitrary and can thus be taken large 
enough to stabilize the fiat directions. 

We next turn to the constraints on the potential at the electroweak scale. In what follows 
we will use just the tree level potential Vo(A). This approximation, discussed in detail in 
Appendix C, is motivated by the fact that we are not concerned with precise predictions 
for the various masses, but rather with the characteristic hierarchies in the spectrum, with 
the rough bounds on the various parameters and with comparing the naturalness of various 
parameter choices. 

First, the scalar configuration with {Hu^) ^ 0, and all the other fields at zero, should be 
a local minimum. This will be the case if we impose that the MSSM parameters, evaluated 
at the electroweak scale, satisfy: 



where as before m| = — 2m^ and m\ = mfj + m|,. (Notice that in the above equation 
we have neglected any contribution to the sfermion masses coming from the Higgs VEV. 
We have also ignored the phenomenological bounds on these masses, which yield somewhat 
stronger constraints: m] > m|. However we stress once more that, for the purpose of studying 
the spectrum hierarchies and the naturalness of different scenarios, the above contraints are 
sufficient. Indeed, in most situations we will end up with sfermion masses well above mz-) 

A second class of constraints is needed to avoid having other minima with electroweak- or 
color-breaking VEVs of order ms- Such minima can arise, even for positive sfermion masses, 
from the destabilizing effect of the trilinear terms in the scalar potential. These are given 
by the soft A-terms and also by the trilinear terms proportional to fi in the supersymmetric 
part of the scalar potential. In what follows we will mainly be concerned with necessary 
constraints, and will not enter into a comprehensive discussion of the sufficient ones. Let 
us consider the effect of A terms first. These were discussed in Ref. |^ where a necessary 
condition to avoid unwanted minima was given: ml + ml + m^ > \A\'^/3 (where a, b, c represent 
any three fields having a Yukawa coupling A, and A is the corresponding soft-SUSY-breaking 
trilinear coupling). When this condition is not satisfied, there is a color- and charge-breaking 
minimum with energy density ~ — /A^. In the case of a light fermion this vacuum is 
considerably deeper than the usual Higgs one. For the top quark, A is sufficiently large, 
and ms is often assumed to be sufficiently small, that this extraneous vacuum is not deeper 
(and typically shallower) than the ordinary vacuum. This is why the A-term requirement 
is usually not applied to the stop. However, when there is a hierarchy ms ^ mz, the 



ml{A ~ lO^GeV) > 



(24) 




(25) 
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extraneous minimum, when it exists, is indeed parametrically deeper than the 0{g'^v'^) Higgs 
minimum, so the necessary condition given above must be apphed also to the soft parameters 
of the third generation. Similar arguments can be made for the trilinear n terms, though we 
are not aware of previous discussions in the literature. Now the triplets of fields in danger of 
developing expectation values are those entering the various /i couplings: {Hu, L, r), (if^, Q, b) 
and [HfiiQit) (where the last member of each triplet is the SU(2)-singlet scalar field). For 
instance, along the direction = (L) = (r) = the potential is given by 

V = {ml + m| + m\m^ - 2^iXM' + \l\ct>\' (26) 

where = + (^^f + (7|)/2. To avoid a minimum away from the origin in the above potential, 
H must not be too bigQ: 

< \l (ml + m| + m?) . (27) 
Notice that, because of the D-term contribution to A^, the bound (27) is irrelevant when 



the Yukawa couplings are small (namely for sfermions of the first two families, or even r 
and h at small tan/5.) When /z is somewhat above this bound, an unwanted minimum with 
V ~ — /i'^A^/A^ is present. Again, for scenarios with a hierarchy ^ ^ m^, this new vacuum is 
much deeper than the correct one. 

In the course of our study we have verified that the positivity constraints of Eqs. (^) and 
( pSD are always stronger than those coming from the trilinear A and ^ terms, at least for the 
parameter ranges of interest to us. Thus, while important in principle, the instabilities arising 
from trilinear terms in the scalar potential do not impose any constraints in practice. 

Next, we examine the evolution of the Lagrangian parameters down to the electroweak 
scale. The form of the RG equations dictates that the soft-breaking scalar mass-squared 
parameters at the electroweak scale are linearly related to their initial values, to the square 
of the GUT-scale gaugino mass M1/2, and to the square of the /i parameter. (In fact there 
are also terms proportional to the GUT-scale values of the A parameters, namely oc A^ and 
oc As shown in Appendix A they can be neglected unless Aq is at least an order of 

magnitude bigger than the other GUT-scale parameters.) Thus the constraints of Eqs. (pSj) 
and Eq. (|2^) , when saturated, constitute a set of eight hyperplanes in the space of initial 
scalar mass-squared parameters. When M1/2 = = 0, the low-energy masses mf are just 
homogeneous linear combinations of the GUT-scale parameters Mf , so the various constraint 
planes determine a cone — or rather, technically, a pyramidal surface — within which those 
constraints are satisfied. Such a cone, drawn in the 3-dimensional space of initial parameters 
for the minimal SO(IO) theory as discussed below, is shown in Fig. 6. For finite gaugino mass, 
A or /X, the constraint planes are shifted by finite amounts. If there was any allowed solid 
angle for M1/2 = A = /i = 0, the new allowed region will be a truncated cone shifted from the 
origin. If there was no allowed solid angle for M1/2 = A = /i = 0, turning these parameters 
on can allow a finite (hyper-)polyhedron. In the absence of running, that is, if the constraints 

^Indeed one can find more general constraints by considering an arbitrary direction in the [Hu, L, fc) space. 
We are not interested here with such a general study — all we want to point out is that fi cannot be much 
larger than the sfermion masses. Similar considerations apply to the A terms. 
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Fig. 6. The allowed "cones" in the space of scalar GUT mass parameters in the minimal 
SO(IO) scenario with exact PQ and R symmetries, within which electroweak symmetry 
is correctly broken. Only the dominant (planar) constraints are shown: m\ > 0, > 
and m? > 0. Note the focusing of the solid angle for increasing Xq, a consequence of the 
exponential homogeneous evolution. 

of Eqs. p5| ) are evaluated at the GUT scale, the cone (or polyhedron) they determine spans 
a solid angle of order unity. As the parameters in Eqs. (|25|) are evolved to lower energies, the 
planes turn about the origin and the cone changes; it may even close completely, in which 
case proper electroweak breaking becomes impossible. We will of course be interested in the 
cone evolved all the way to the electroweak scale; it is useful to remember that it is drawn in 
the space of parameters in the effective GUT-scale Lagrangian, and its boundaries correspond 
to those GUT-scale parameter values which lead to the vanishing of particular scalar masses 
at the electroweak scale or of ml at 10^ GeV — in other words, it encompasses the GUT-scale 
parameters which would lead to proper breaking at the electroweak scale. A narrow cone 
means that it is difficult to find GUT-scale parameters which will lead to a low-energy world 
similar to ours. 
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6.2 The Homogeneous Evolution 

To understand the evolution of the cone, for the time being we restrict our attention to the 
homogeneous parts of the differential RG equations for the soft parameters, neglecting the 
gaugino, A and /i contributions. This not only illuminates the functional behavior of the 
solutions, but is also directly relevant for the case we have so far espoused, the PQ- and 
R-symmetric one in which the /i and A parameters and the gauginos are much lighter than 
all the other masses (except m^). In this maximally-symmetric case the evolution is driven 
by the Yukawa couplings, which in turn depend on Xq (and pa)' 

• For small Xq there is little evolution, and the cone remains wide. Here there are no 
generic difficulties in satisfying the constraints for a wide range of initial parameters. 
Whether such a range is available in particular GUT models is a question that will be 
answered in the next section. 

• To understand the behavior for large Xq, it is useful to change basis, considering certain 
fixed linear combinations A1m| (where is a constant matrix) of the seven soft- 
breaking parameters. Appendix A contains the solutions, in this new basis, of the 
RG equations for the seven soft scalar masses, as well as for the three A terms and 
the B parameter, in terms of integrals over only the gauge and Yukawa couplings. The 
matrix Ai is chosen to separate the homogeneous part of the seven scalar mass RG 
equations into two classes: three of the linear combinations, denoted Xt, Xf, and X^, 
evolve essentially multiplicatively, contracting exponentially as they evolve down to the 
electroweak scale. The other four, /i,2,3,45 are approximately invariant. (Actually one 
linear combination of the Jj evolves slightly due to the hypercharge D-term S discussed 
in Appendix A, while three other independent combinations are truly invariant; but for 
the present purposes we can neglect S and consider all four Jj to be invariant.) Thus 
we find that, for large Xa, the three Xi are exponentially suppressed at the electroweak 
scale relative to their GUT-scale values, and hence, gencrically, also relative to the four 
invariants /j. In the limit of very large Xg, at which the Yukawa couplings tend to their 
"fixed points" at low energies, the constraint equations constitute seven constraints in 
only four variables Jj: 

< m,^ = M-i ^ ) ^ M~' ^~^^'^''^G ^r.M-'(^fy (28) 

Given the value of Ai in Appendix A, it is simple to check that, in the limit A — > oo, 
there is only a trivial solution: Jj = 0. This solution is highly nontrivial, however, in 
the sense that it requires adjusting the soft masses at the GUT scale such that they 
just cancel when combined into the 7^. At finite A the degree of cancellation needed is 
just the degree to which the X^ are suppressed; schematically, then, we must adjust the 
GUT-scale parameters to satisfy 

\h\ < e-^ I Xg, 2 = 1,2,3,4. (29) 
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How much tuning is actually required in order to satisfy Eqs. (^5D ? We will determine 
numerically the allowed region of parameter space under various GUT-scale assump- 
tions. Generically, one would expect a factor of e"^-^'^ '^^ for each invariant which 

needs a cancellation, hence ( e ^ J 1 . However, since the leptonic Yukawa coupling 

is considerably smaller than the other two over much of the RG evolution, and since 
one of the invariants and one of the X's are essentially determined by leptonic masses, a 

better estimate is ^e"^-^'^ '^^^ . To alleviate this tuning problem, we should stay away 

from the fixed-point regime of large Yukawa couplings, or relax the PQ or R symmetries 
(but then e^e^ is made worse). Away from the fixed point, a splitting Xf > also 
helps. In the end, of course, the question itself only makes sense in a particular context: 
so far we have assumed there is no theoretical bias towards any relationship between 
the initial parameters, except for the approximate PQ and R symmetries. In a GUT 
context there will be some definite biases. 

The large-Ac fine tuning for the symmetric case is illustrated graphically in Fig. 6, assum- 
ing the SO(IO) scenario discussed below in which only one contracted direction X and two 
invariants /i,2 play a role. The RG equations for the seven soft masses are solved for three 
cases: Xq = (that is, without running), an intermediate value of Ac, and a large Xg- The 
figure shows the cone in which the electroweak-breaking constraints of Eqs. (p5D are satisfied, 
drawn in the space spanned by the three linear combinations of GUT parameters Xq, h and 
l2- Note that only a small number — three or four — of the constraints are the decisive ones, 
and once they are satisfied all others are as well; in this case, they are m| > 0, m\ > and 
m? > 0. We see that near the fixed point, the cone closes up0 around the ray Ii = I2 = 0, 
meaning that if we don't tune the GUT parameters to lie in this ray, then the soft-breaking 
parameters at the electroweak scale will not satisfy the constraints. Then either electroweak 
symmetry will not break (m| < 0), or the potential will not be bounded from below {m\ < 0, 
leading to dimensional transmutation at a scale much larger than the SUSY scale, and also 
to tan/5 ~ 1), or an electromagnetically-charged scalar will acquire a VEV. 

The contractions may also be seen analytically in terms of sum rules, which are particular 
linear combinations of the electroweak-scale soft-breaking masses having only positive coef- 
ficients and chosen to be independent of all the invariants. Linear combinations of the Xj 
with positive coefficients can give such sum rules, for example Xt -t-X^. These sum rules have 
the property that on one hand they are phenomenologically constrained to be positive (since 
mf > 0), but on the other hand they are driven to zero by the RG equations as Xg increases. 
Schematically, we have for our example 

< 2m| + m? + m? + = Xt + Xb ~ e"^ / (X^ + Xb)^ . (30) 

Since each of the mass terms in this equation should be positive, for large Xg they must each 
be made to evolve towards zero at low energies. Now, each mass term can itself be expressed 

^Actually, as evident from the figure, the cones close up before they reach the ray Ii = I2 = 0. This 
"premature focusing" is a property of the specific GUT boundary conditions and wih be discussed below. 
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as a sum of the Xi and the If, since the Xi are exponentially reduced while the li remain 
invariant, the soft masses can only evolve towards zero if the various /j are tuned to be small 
already at the GUT scale. In fact, the sum rules embody the same information as the planes 
in the constraint equations (|25| ) and in Fig. 6. In particular, there is a dominant sum rule 
corresponding to the innermost set of planes which define the constraining cone; for the case of 
Fig. 6, this dominant sum rule is m| + ^{m\ + m?). It indicates which are the masses closest 
to saturating the constraints, and hence which are typically the lightest. Notice also that 
all the soft scalar masses except the sleptons appear in the sum rule of Eq. (0), indicating 
that essentially all these parameters contract for large Xq- The slepton masses also contract, 
according to another sum rule, but to a lesser degree. 

A very useful graphical way to describe the allowed domain in the space of initial pa- 
rameters is to project the constraint planes onto the (hyper-)plane spanned by the various 
invariants. For the SO(IO) scenario illustrated in Fig. 6 and described in detail below, the 
result is a set of lines in the plane of {Ii/Xg, h/Xc) ■ To normalize the axes, we define a 
"typical soft scalar mass" Mq = ^Xq = ^M^q^ + |M^gg, and use the more direct Lagrangian 
parameter Mx (see below) instead of l2- The horizontal and vertical axes are thus shown in 
units of M^/Mq and /i/Mq, respectively. Each line forms the boundary of the half=plane 
where one of the is positive. Fig. 7 shows all eight lines, and emphasizes the region allowed 
by the constraint equations (p^ and (^51) , for three values of Xg, and for various values of 
Xf /Xf. In general the allowed region is a polygon; for the PQ- and R-symmetric SO(IO) case, 
it is usually a triangle bounded by the lines corresponding to the three masses which appear 
in the dominant sum rule m| + ^{m^ + m?). Near an edge of the triangle, the corresponding 
mass parameter is much smaller than all the others. The hatched region in Fig. 7(a) is where 
m| is ~ tan/3 ~ 50 times less than m\, which is the favored scenario (as discussed in more 
detail below). Similar regions are indicated for the other triangles. 

6.3 Evolution and Natural Selection 

From the previous subsection, we conclude that if /i and are chosen much smaller than the 
typical soft scalar masses so PQ and R are approximately valid, then we expect the allowed 
triangular area in the space of GUT-scale scalar masses to be small if ^ 0.7 and = A^ 
(the focused case), and large otherwise. In fact, as we show below, both SO(IO)- and SU(5)- 
type boundary conditions on the scalar masses result in premature focusing: the triangles 
close up for finite Xq values, not far above unity. Within the allowed triangle, therefore, a few 
particles — namely those which bound the triangle itself — are very much lighter than the rest. 
(In contrast, for Ac — > cxd focusing all masses must become very light.) 

Now, if there is not much focusing, then all the scalar masses are comparable throughout 
most of the triangle, while in a narrow (~ 1/ tan (3) strip within that triangle m| is ~ 1/ tan/3 
times lighter than m\. Tuning the GUT scalar parameters to lie within this strip suppresses 
the large b ^ s'-f amplitudes, and allows us to generate large tan (3 with no further tuning of 
B and without violating LEP bounds by taking /i^ ~ ^1/2 ~ "^a/ there is significant 

focusing, then within the small allowed triangle the particles which bound the triangle have 
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= 0.4 = 0.7 = 1 




Fig. 7. The allowed regions for the same scenario as Fig. 6, but projected into two 
dimensions by dividing by the typical GUT squark mass Mq. The GUT-scale Yukawa 
coupling is increased from left to right, and both exact Yukawa unification (black areas, 
solid lines) and approximate unification (light gray areas, dotted lines) are considered. 
For Ag = 1 we also show the perturbed case X^/Xf = 1-1 as the dark gray area and 
dashed lines. The lines delineate the half plane in which the corresponding mass-squared 
is positive and hence acceptable. In the diagonally-hatched regions <C rn\, as 
discussed in the text. Finally, the vertically-hatched triangles are the allowed areas 
assuming exact unification and a right-handed neutrino at Mn = 10^^ GeV. 

masses — indicated collectively with mg — which are much lighter than the others, so a hier- 
archy ml ]^ -C m^jj is invariably present. As we will see in Tables 3 and 4, m| and are 
always among the light masses, while m~ is always one of the heavy masses. Hence throughout 
the triangle the dangerous b ^ sj amplitudes are somewhat suppressed, depending on the 
degree of focusing. There is now no naturalness criterion to distinguish between the following 
situations: either the GUT scalar parameters are tuned to lie within a narrow strip in this 
triangle, resulting in m| ~ ^'a/ tan/5 ~ '^ol ^"^^ requiring no further tuning of B since we 
can select /x^ ~ ^1/2 ~ "^1 ~ "^a/ tan/? and meet the LEP constraints; or the GUT scalar 
parameters are not adjusted to be in the strip, so m| ~ m\ ~ ml j^, but B must be tuned to 
~ 1/ tan/? of its natural value because LEP requires us to select /x^ ~ ^1/2 ~ "^1 ~ "^a- To 
summarize: when the PQ and R symmetries hold at the GUT scale, if there is no focusing 
then they hold at the electroweak scale and the overall tuning need only be ~ 1/ tan/3; but if 
there is strong focusing then they can either hold or not hold at low scales, and in either case 
the need for tuning is greater than 1/ tan/3 by the degree of focusing. 
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Strictly speaking, the exactly PQ- or R-symmetric scenarios are never allowed by LEP 
limits, so we should in principle always take fi and M1/2 into account in evolving the cone to 
low energies. But this would not qualitatively change the discussion. Clearly, approximating 
fj, or Ml/2 by zero is valid when there is no focusing and the symmetries are approximately 
valid at all scales. But even when there is focusing and some of the particles end up with 
small masses mo,L, introducing finite /i or M1/2 comparable to those masses will only change 
the allowed triangle area and the light masses by 0{1)] in fact, if /i or M1/2 are then increased 
beyond mo,L, the focusing is alleviated such that the new value of mo,L is again comparable 
to /J, or Ml/2. 

If the PQ or R symmetries are not approximately valid at the GUT scale, the planes m? = 
which delimit the allowed volume are shifted, and the focusing is alleviated. (The exception 
is when fi or M1/2 are too big in certain GUT scenarios — then the triangles once again close 
for finite values of Xq-) Irrespective of focusing, the tuning is worse than ~ l/tan/3, as 
demonstrated schematically in Table 2. However, we will reserve judgment on the breaking of 
PQ and R symmetries until after we have studied the SO(IO) and SU(5) symmetry constraints 
in the scalar mass sector. 



7 Grand Unified Soft Masses 

7.1 SO(IO) and SU(5) Boundary Conditions 

We now address the question of what values of initial parameters arise from theories at the 
GUT scale. Since the idea of Yukawa unification is largely based on symmetry principles, 
it behooves us to consider the implications of those same symmetries for the various soft 
SUSY-breaking parameters. We have already employed, in our solutions to the RG equations, 
the assumption that gauge coupling unification is accompanied by gaugino mass unification, 
in order to reduce the number of independent gaugino mass parameters to one. Let us now 
examine the implications of SO(IO) or SU(5) symmetries for the soft-breaking masses. As a 
special case, we comment briefly on the universal scenario. We also examine the threshold 
corrections due to a light right-handed neutrino. 

Consider first the simplest SO(IO) scenario, in which both light Higgs doublets originate 
from a single IOh multiplet, or more generally any SO(IO) model in which all the GUT fields 
from which the light doublets arise have degenerate soft masses. When SO (10) breaks, in 
general both D and F terms could split the scalar masses in a single SO(IO) multiplet. D-term 
splittings are generically present because the rank of the gauge group is reduced, but 



F-term splittings are more model dependent and need not arise. For example, when SUSY 
breaking is communicated from a hidden sector only via gravitational interactions, the soft 
terms are very constrained EOl. This property leads, in the minimal missing- VEV models 



4l[| , to exclusively D-term splittings. Therefore when we refer to SO(10)-type boundary 
conditions we will only include D-term splittings, whereas more general F-term splittings 
will be encompassed by the discussion of SU(5) boundary conditions, or when necessary by 
the completely general discussion. Accordingly, in SO(10)-type models the seven soft-breaking 
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masses are determined at the GUT scale by only three soft-breaking parameters: the soft Higgs 
mass Mio^, the third-generation soft squark and slepton mass M163, ^ mass parameter 
Mx from a D-term that is left over when the U(l)x symmetry in SO(IO) is spontaneously 
broken. Recall that the rank of SO(IO) is higher by one than that of SU(5) and of SU(3) ® 
SU(2) (g) U(1)y, thus SO(IO) D SU(5) (g) U(l)x, where the generator of U(l)x is proportional 
to 3(B — L) + 4T3R. One common way to break this U(l)x and reduce the rank to that of the 
standard model is to introduce a pair of 16-dimensional Higgs representations having GUT- 
scale masses, 16h and 16^, which acquire VEVs in their 'V/j" components, thus preserving 
the SU(5) symmetry. When U(l)x breaks this way, its D-term acquires a VEV proportional 
to the difference of the soft masses of the 16h and 16h- This D-term then contributes to 
the soft masses of the fields which couple to U(l)x in proportion to their U(l)x charges. To 
quantify this contribution, we define a mass parameter 

M'x = Ml,^ - M^^ . (31) 

(Note that this definition differs by a factor of 10 from the definition in our previous work 
|1T6[| . We chose the present definition because it is closer to the fundamental parameters of the 
grand-unified theory, and so is on the same footing as Mio^ and M163.) presence of 

Mx which allows the up- and down-type Higgs masses to be split at the GUT scale in almost 
any SO(10)-unified scenario, and thus greatly facilitates proper electroweak breaking. 

In SU(5), of course, the two Higgs doublet superfields originate from different representa- 
tions, 5h and 5^, so their soft SUSY-breaking masses are generically split. This is also the 
case in some SO(IO) models, for example when the light doublets are mixtures of different 
GUT multiplets having different soft masses, or even when they both lie in the same SO(IO) 
multiplet but the soft terms are general enough to induce F-type splittings. Note that in this 
second case the Yukawa couplings are indeed exactly unified. For brevity, however, we will 
call any such boundary conditions on scalar masses SU(5)-type. The seven soft masses in 
SU(5)-type models are determined by four parameters: the two Higgs masses M^^ and Mg^, 
and the soft masses for the two representations which contain the third-generation MSSM 
squarks and sleptons, M^^ and M103. For comparison with the SO(IO) case, we will recom- 
bine these four parameters into the same three combinations which occur in SO(IO) plus an 
extra degree of freedom Msu(5), as follows: M^^^ = i(M|^ + M^J, Mi^.^ = fM^^^ + iM|, 
M'x ^ liMl - M|J + liMi,^ - Ml), and M^^^,^ ^ \{Ml^ - M^) + \{Ml,,^ - M^). (To 
reiterate: in an SU(5) context, the four quantities Mio^, M163, and Msu(5) should be 
regarded just as convenient linear combinations of the underlying soft masses.) 

With these redefinitions we may write the seven soft masses at the GUT scale in either 
SO(IO) or SU(5) as follows: 

Mh = Ml^^ _lMl + iMs%5) (32) 

Ml = Ml,^ +1m^-1Ms%5) (33) 

1 .,0 1 



Klf = ^163 + + 4^lu(5) (34) 
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= M^^ea - - ^M|u(,) , (35) 

(Capital letters denote the parameter values at the GUT scale.) Their linear combinations Xi 
and then take the form: 



X? = 








+ ^lu(5) 


(36) 


X? = 




+ 2M2,3 




~ ^SU(5) 


(37) 


X? = 




+ 2M2,3 




~ ^SU(5) 


(38) 


h = 


2 Ml^^ 






+ 2M|u(5) 


(39) 


h = 






5 ^ 


+ ^^lu(5) 


(40) 


h = 




Ml,, - 


-Ml 
10 ^ 


+ i^lu(5) 


(41) 


h = 








~ ^SU(5) ■ 


(42) 



A common assumption in much of the previous work on unified supersymmetric theories, 
and in particular in large tan/3 models is a "universality" of soft SUSY-breaking scalar 
masses at the GUT scale. Under the universality assumption, Mio^ = Miea while Mx = 
-^su(5) = 0. We have not made this assumption because we do not expect it to hold at the 
GUT scale (whether or not it is a good approximation at the Planck or string scales) and 
because it requires |l^ tremendous fine-tuning of parameters to achieve proper electroweak 
symmetry-breaking with large tan (5. We will have more to say about this case in the discussion 
in Sec. 

A much better motivated assumption is that the third-generation right handed neutrino, 
that is, the electroweak-singlet superfield which couples to vl through a Yukawa coupling 
with Hui has a Majorana mass Mn smaller than Mqut- Theoretically, such a mass may 
arise^ from a high-dimensional operator induced at some scale M' > Mqut (such as the string 
or Planck scale) and therefore be suppressed by a power of 1/M'. Phenomeno logically, a lower 
Mjv leads through the see-saw mechanism to a higher mass for the observed neutrinos, which 
may then explain various cosmological and astrophysical puzzles. In any case, as long as M^r 
is not too far below Mgut, its effects can be approximated by threshold corrections to the 
Yukawa couplings and to the soft SUSY-breaking scalar masses. If the SO(10)-type boundary 
conditions Xf = X% = Xg are valid, we can calculate the size of the corrections: 



2 ^ 87r2 



^This happens, for instance, in the absence of 126 + 126 SO(IO) Higgs multiplets and of R-odd gauge 
singlets. Note that in order to preserve the vaHdity of perturbation theory it is better to avoid the 126 or 
bigger representations, since they contribute a large positive term to the gauge /3-function [M. 
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(43) 



The consequences will be discussed in the following section. 



7.2 SO(10)-type GUT Masses 



We first consider the symmetric, minimal SO (10) scenario: (PQ- and R-) symmetric in that 
H, A, B and the gaugino masses are negligible relative to the various squark, slepton and 
pseudoscalar masses; minimal in that the two light Higgs doublets are contained in a 10^^ of 
SO(IO), so = (up to small threshold corrections); and SO(IO) in that the soft-breaking 
masses are those that arise in a minimal SO(IO) theory, hence Msu(5) = 0. Three independent 
dimensionful GUT-scale initial parameters specify the electroweak-scale consequences, and we 
choose them (as noted above) to be (= \MIq^ + fM^gs = = l^b = = 1^^)^ 
Mx {— 5/2) and 2MIq^ — 3M|gg = 7i. The third invariant is a linear combination of Mq and 
the first two invariants; it is not very constrained, since it is mostly associated with Xr, which 
does not contract much. The fourth and last invariant vanishes identically by the SO(IO) 
symmetry. The evolution from the GUT scale to the electroweak scale, and the contraction 
of the allowed triangle towards small Mx and Ji as A^ increases, were illustrated in Fig. 6, 
or in projected form (in the space of Mx/Mq and Ii/Mq) in Fig. 7. Three features are worth 



1. In the general case of seven independent initial soft mass parameters, we expect that 
the allowed region closes asymptotically as Ac ^ 00, when all four invariants must be 
set to zero to allow the three Xi to contract. Graphically, this means that the eight 
planes (or lines) corresponding to mf{A = rriz) = and m|(A = 10^ GcV) = all 
cross at one hyper-ray /i,2,3,4 = (namely the origin in the projected space). Instead, 
we see that in the symmetric minimal SO (10) case the allowed region of parameter 
space is bounded by the planes (or lines) m| = 0, m\ — and — which cross 
prematurely, closing the allowed triangle for a finite value of Ag- The reason is that 
SO (10) boundary conditions are compatible with /i,2,3,4 = only if all initial masses 
vanish, a trivial and uninteresting scenario. For nonvanishing Mq. the restrictive SO(IO) 
boundary conditions can only be satisfied for sufficiently small Xq- Just how small? 
Consider the sum-rule- like combination m| -|- |(m^ -|- m~), which was chosen so that 
the invariant part is a negative definite quantity: (— 2Ji — 7/2 — 14/3 -|- 1474)/183 = 
~A/61Ml The contracting part {-97Xt + 216X5 + 7X^)/183 starts out positive, with 
a value 126/61Mo, and contracts monotonically to zero as Xq increases. Hence the sum 
^1 + K^i + ^l) = (-97Xt + 216Xb + 7X^)/183 - 4/61M2 must vanish for a finite 
Xg- The critical value turns out to be Ag — 1-2. Therefore in this minimal symmetric 
SO(IO) case it is important that Ag is not only small, but in particular is well below 
~ 1.2, for there to be a significant allowed region of parameter space. [If we wish to 



noting: 
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include the small contribution to this "sum rule" of the hypercharge D-term, we can 
again find the proper combination of masses in which the invariant part is negative and 
proportional to Mq . The result is: 

but there is hardly any change in the conclusions.] 

2. When the parameter space contracts, it does so around a nonzero value of M'^/Mq. This 
highlights the important role of the U(l)x D-term mass parameter in allowing both a 
negative mfj and a sufficiently positive m|,, even though the RG equations drive m|) 
down more than m^. 

3. Notice that the allowed cone in the PQ- and R-symmetric SO(IO) scenario of Fig. 6 
is confined to the Mq > half plane. This result was used when dividing by Mq to 
determine the allowed area in the projected space of Fig. 7. We have found that, in 
order to satisfy the low-energy constraints of Eqs. (^) and Eq. (p4D, Mq must always 
be positive for all values of Xg and \f /X^ under consideration, and for SU(5) as well as 
SO(IO) boundary conditions on the soft scalar masses, as long as the PQ and R sym- 
metries hold. To meet the low-energy constraints with negative Mq inevitably requires 
very large values of I/i^/MqI and sometimes of |Mfy2/M?h since these would require 
very delicate fine-tuning of the GUT-scale parameters, we will assume Mq > for the 
remainder of this paper. 

From Fig. 7 we can infer some properties of the universality assumption when PQ and R 
are vahd. The asterisk (*) indicates the coordinate (M|,/i) = (0, — Mq ) corresponding to 
universal scalar mass boundary conditions at the GUT scale. Proper electroweak symmetry- 
breaking occurs only if and are widely split, and then only for intermediate values 
of Ac. Furthermore, to meet LEP constraints in the approximately PQ- and R-symmetric 
scenario, the value of Xg must be tuned to achieve m| ~ /^"^^ ^ ^o- If on particular we 
set /i^ ~ m~, ~ m\/ tan/3 and tune Xg with a precision ~ 1/ tan/3 to get m| light enough, 
we end up with the minimally fine-tuned scenario at the lower- left corner of Table 2. Thus 
the PQ- and R-symmetric universal case is allowed and only minimally tuned (via Xg) if the 
Yukawas are widely split at the GUT scale. 

Fig. 7 also illustrates the effects of threshold corrections due to a right-handed neutrino 
with a mass M^ -C Mqut- The vertically-hatched triangles in Fig. 7 shows the area allowed 
when Mjv ~ 10^^ GeV. The ~ 5% correction to the Yukawa couplings lowers Xf and so would 
reduce the allowed area, but the ~ 40% negative correction to M^ is far more significant and 
increases the area. The result is that the area remains ~ 0.1 even at Xg — 1.0, comparable 
to that of the dark-shaded triangle [shown only in Fig. 7(c)] which would arise from a typical 
~ 10% threshhold correction to Xf/Xf. Notice, however, that the neutrino effect, in addition 
to being well motivated, also favors a small Mx ~ 0; in other words, the AM^ shift can 
substitute for the shift produced by the U(l)x D-term which was needed for proper radiative 



gQG ln(MGUT/mz) 

1 + Sac In(MGUTMz) 



(44) 



31 



symmetry breaking. A vanishing could conceivably be achieved naturally by means of a 
symmetry. Finally, note that a light right-handed neutrino threshold is not sufficient to allow 
for the PQ- and R-symmetric universal case when = A^. 

In Table 3 we display some characteristics of the SO (10) scenario for various values of 
Ag and of Xf /Xf. The heavy-dashed boxes correspond to the PQ- and R-symmetric case: 
the three on the left are for the minimal Higgs choice (A^ = A^^), while the three on the 
right allow for large Higgs mixings at the GUT scale (A^ = 2Xf). The top entry in each 
box gives the area of the allowed triangle using the coordinates of Fig. 7, namely M^/Mq 
and Ii/Mq. Also shown in some interesting cases are the larger areas that would result 
from a slight Yukawa sphtting (Af = ^-^^b) some shght mixing or threshold effect. 

Note that the area decreases rapidly as Ag increases, indicating the aggravated need for fine- 
tuning of the GUT-scale parameters. Of course, the value of the area depends on the choice 
of coordinates and the metric, which are to some extent a matter of taste. We use these 
particular coordinates because we expect them to be a priori of order unity, and so if the 
triangle area is much smaller than 1 then some tuning is apparently necessary. A crossed-out 
box indicates that the corresponding parameter choice leads to a value for Sj, incompatible 
with bottom-tau unification and the low-energy values for and (though some of those 
boxes are nevertheless ffiled in for reference). 

The full superspectrum is completely determined, up to an overall scale, by choosing a 
specific point within the allowed triangle. For the PQ- and R-symmetric scenarios we are 
now considering, the middle of each box indicates the superspectrum that is typical near the 
light-m| portion of the triangle, shown as the hatched region in Fig. 7. [Recall that, while 
this region is ~ 50 times smaller than the triangle and hence by definition requires that much 
fine-tuning, it leads to a hierarchical spectrum which is both phenomenologically allowed and 
requires no further tuning to achieve acceptable tan/3 and V{h — > 57). Furthermore, far away 
from this region, niz is greater than fi and and hence is in confiict with LEP.] In most 
cases all the superpartners and the pseudoscalar Higgs have similar masses. However, when 
Ag is large and Xf ~ A^, the allowed triangle is small and therefore its bounding particles, 
the pseudoscalar Higgs and the SU(2)-singlet bottom squark, are somewhat fighter than the 
other particles. And when Ag is large but there is a large top-bottom Yukawa splitting, the 
SU(2)-singlet stop becomes relatively light. The last item in each box is the list of masses 
which can vanish simultaneously with m^: for the symmetric SO(IO) case, they are always 
mz, rriA and m^, as discussed above. In other words, we may choose parameters at the corner 
of the allowed triangle such that tua and mz are much lower than all the other superpartner 
masses, or and mz are much lower than the others. 

Fig. 8(a) shows contours of fixed allowed area as functions of the size of the GUT-scale 
Yukawa coupling Ag and the amount of top-bottom splitting splitting Xf /X'^. The sharp 
bends occur when the m? > constraint becomes more restrictive that the m? > constraint, 
SO the rate at which the triangle closes is determined by the evolution of m|, m\ and m| 
rather than of m|, m\ and m~. Note the dramatic decrease in area as the maximal value of 
Ag is reached for fixed Xf/Xf — ^this is the premature focusing implied by SO(IO) [and SU(5), 
as we shall see] boundary conditions on the scalar masses. The large-AG, small-splitting region 
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Table 3. The characteristics of exact and approximate Yukawa unification with SO(IO)- 
type boundary conditions on the scalar masses. For each of the six choices of Xf /X^ 
and Xq, four boxes are shown, corresponding to the presence of absence of approximate 
PQ and R symmetries (with fi/Mo and M1/2/MQ shown respectively in square brackets). 
The first entry in each box is the allowed area, the second is a typical spectrum, and the 
third lists the masses than can be decreased simultaneously with niz- See the text for 
further details. 



clearly requires very precise adjustment of the GUT-scale parameters. 

We next relax either the PQ or the R symmetry, or both, and once again ask for the 
allowed region in the parameter space of soft scalar masses. For fixed /i and M1/2, the planes 
which delimit the allowed region are now shifted by fixed amounts ~ /i^ and ~ M^^, so 
that they no longer intersect at the origin. We are actually most interested in the relative 
quantities fi/Mo and M1/2/M0, that is, the amount of PQ and R breaking relative to the 
other (soft scalar) mass parameters. So for fixed n and M1/2, we should consider various slices 
of constant Mq in the scalar mass parameter space. (Such slices are in fact the projections 
shown in Fig. 7.) Small PQ and R breaking corresponds to looking at large- Mq slices: at 
such large distances from the origin, the small displacements of the planes are insignificant, 
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Fig. 8. Contours of constant allowed areas or volumes for SO(IO)- or SU(5)-type bound- 
ary conditions, respectively, assuming exact PQ and R symmetries. Notice that small 
Ac and large Xf /Xf values are favored, since they alleviate the premature focusing of 
the homogenous RG evolution. 



and the allowed portion of the slice is essentially the same as in the symmetric case. Larger 
PQ and R breaking correspond to slices of smaller Mq, nearer the origin: the allowed portion 
of the slice may be large relative to the distance from the origin, in which case there is no 
fine-tuning, or the slice may not even intersect the allowed region, in which case those values 
of /i/Mo and Mi/2/Mq are not allowed. 

We must also ensure that, within the allowed region, Sb is within the range allowed by 
bottom-tau unification. [To calculate 6b for a particular choice of /i/Mo and M1/2/M0 and a 
given value of as(m^), we first determine the range of values of m| + m~ within the allowed 

region of parameter space, and then use Eq. if any resulting 6b is acceptable, we allow that 
choice of /i/Mq and M1/2/M0.] In Fig. 9 we show the values of /x/Mq and M1/2/M0 which lead 
to proper electroweak symmetry breaking (i.e. there is a finite allowed triangle) and acceptable 
bottom-tau unification, for three representative values of Ac, two values of a<j(m^), and either 
\f = (shaded) or \f = 2A^ (hatched). Recall that M1/2 is the gaugino mass at the GUT 
scale (which happens to roughly equal the wino mass: = g^/ gQMii2 — O.85M1/2) and Mq 
is the typical scalar mass also at the GUT scale, while /z is evaluated at the electroweak scale. 
The allowed regions are all roughly "L" -shaped. At their top and on their far right (when 
applicable), they are cut off by the requirement of proper electroweak symmetry breaking, 
while on their lower-left and upper-right sides they are bounded by the limits on 6b- 
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Mi/2/Mo as = -115 




as = .125 




Fig. 9. The regions of GUT-scale PQ and R breaking allowed by the constraints of 
proper electroweak breaking and bottom-tau unification. The gray areas give the ranges 
of ^/Mq (where /i is actually evaluated at a scale mz) and M112/MQ assuming exact 
Yukawa unification, while the hatched regions assume A^/A^ = 2. All assume SO(IO)- 
type scalar mass boundary conditions. 



Particular values of fi/Mo and M1/2/M0 are examined in more detail in Table 3. Once 
again we consider three different values of Xq and either = A^^ or Xf = 2Xf, which lead 
to the six major boxes of the table. Each is divided into four sub-boxes: the lower-left one 
is the symmetric case described above; in the lower- right one we break PQ, in the upper- left 
one we break R, and in the upper-right sub-box we break both the PQ and the R symmetries. 
(Note that these correspond spatially to the four corners of Fig. 9, as well as to Table 2.) For 
each sub-box we have chosen an appropriate pair of [yu/Mo, M1/2/M0] values from the allowed 
region of Fig. 9, and have indicated the allowed area of the triangle for those values. When 
choosing these values we avoided the boundaries of the allowed regions, because there the area 
is typically very small and hence (in some sense) unlikely. Clearly these choices are somewhat 
arbitrary, and there can correspondingly be some variation in the spectrum. Notice that when 
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the allowed area is small in the symmetric case, raising /i or M1/2 has the expected effect of 
increasing the area (but recall the price one pays in e^e^), since the focusing described above 
is alleviated. When either the PQ or R symmetries are approximately valid, it is still necessary 
to focus on the light-m^ part of the triangle (see Table 2); when both symmetries are broken, 
much of the triangle is allowed by LEP and no part is selected by naturalness criteria, so 
we have arbitrarily chosen to look at its center. The middle entry in each sub-box indicates 
as before the typical superspectrum for the appropriate region of the triangle. If the various 
squark, slepton and pseudoscalar masses for a given sub-box are all within roughly a factor 
of two, we characterize them as "all comparable", and otherwise we indicate which ones are 
significantly lighter; masses in parenthese are only marginally lighter (i.e. somewhat less than 
half the heavier masses). The bottom entry of each box shows once again the particles which 
can become light simultaneously with the Z. If only one particle is shown, the reason is that 
at the other vertex of the triangle the flat-direction mass is negative at scales between 10^ 
and 10'^ GeV, and hence that vertex would lead to improper electroweak symmetry breaking. 

Returning for a moment to the universal scenario, we recall that if = then the 
universal assumption is incompatible with even approximate PQ and R symmetries. Indeed, 
when the Yukawa couplings are exactly unified the universal case requires a much bigger 



tuning |T6[. The reason is that when the scalar masses are universal the only sources of 
custodial breaking available for splitting mfj from mj^ are the small effects of hypercharge 
and A,-. Moreover, in order to obtain mfj < rn^) the gauginos must be very heavy. This can be 
represented by the low energy relation rn^y — rnfj = ecMy2y where ec is small positive coefficient 
representing the custodial breaking induced by hypercharge. Using this relation, we see that 
proper symmetry breaking, i.e. mfj < < mj^, requires tuning some parameter (Ac, ^"^ or 
My2) to a precision of ec- Moreover, since m\ = ml; + m|, < ecMy2, and since /i ~ C(Mi/2) 
to make m\ > when the gauginos are heavy, we must tune B with a precision ~ 
(/im^/m^)(l/ tan/?) ~ ec/ tan/?. Thus the overall tuning is at least ~ e^/ tan/5 ~ 1/ tan^/3, 
using the rough numerical approximation ~ 1/ tan/5. Some more tuning is required to 
achieve an acceptably small rate for b 57. And finally, with the large 6b corrections that 
result from such a spectrum, the top mass is quite light and is therefore in conflict with the 
recent data on [p6 . 



7.3 SU(5)-type GUT Masses 

We can repeat the above analysis for SU(5)-type boundary conditions, that is, when the soft 
scalar masses need only be SU(5)-symmetric at the GUT scale. There are now four rather 
than three independent initial masses, and we will choose them to be Mq [= |M^q^ + ^MIq^ = 
i(M|^ + + Ml) + iM^oJ, Mi f (M^ - M|J + l{Ml,^ - Mj)], 2Ml,^ - SM^,^ 

[= Mi, + - fM2„3 - |M|J, and Mi^^,^ i(M|^ - M|J + i(M^03 - ^1)]. The 
allowed region in the 3-dimensional projected space of initial scalar mass parameters is now 
a volume bounded by planes, which in many cases is a tetrahedron [corresponding to the 
allowed triangle in SO(IO)]. When the PQ and R symmetries hold, we find as for SO(IO) 
boundary conditions that the asymptotically focused case Ac — ^ 00 cannot be reached, since 
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the four equations /i, 2,3,4 = have no nontrivial solutions. Hence the allowed tetrahedron 
closes for a finite Xq — 1-7, which can be seen by considering the evolution of the combination 
ml + §{mA + ml) + fm? = {-UXt + 39Xb + 15X^)/34 - 6/347^0^ (neglecting the small 
hypercharge D-term). 

Table 4 summarizes the consequences of imposing SU(5)-type boundary conditions on the 
soft mass parameters for various choices of Xq and Xf/Xf, and for unbroken or broken PQ 
and R symmetries, in analogy with Table 3. This time it is the allowed volume, rather than 
the area, which is shown on the first line of each sub-box; also shown arc the sampled values 
of [/i/Mo,Mi/2/Mo]. This volume is also plotted, as a function of Xg and Xf / X^ , in Fig. 8(b). 
The choice of coordinates in the 3-dimensional initial parameter space was as usual a matter 
of taste, so there is no objective way of comparing the allowed volumes in SU(5) with the 
allowed areas of SO(IO). Qualitatively, however, it seems clear from Table 4 and Fig. 8 that 
SU(5)-type boundary conditions require less tuning, mainly because the additional degree of 
freedom Ml^^g-j makes the Higgs splitting independent of the squark masses and hence more 
easily allows mij < mjj without lowering m?~. The middle entry in each sub-box describes 

the typical spectrum near the m| = face of the allowed volume, except for cases of broken 
PQ and R in which the spectrum is shown for a generic central point in the tetrahedron. The 
bottom entry shows which masses are allowed to vanish simultaneously with mz', in general 
there are corners of the allowed region in which two of these may vanish along with mz, but 
which pairs may do so varies from case to case. Finally, notice that in some cases, such as 
Ag = 1 and M1/2 = 0, the SU(5) entry shows that only ttia and nif may be light along with 
mz, whereas the SO (10) entry indicates that only and m^ may do so. The reason is that 
SU(5) allows the larger value /j/Mq = 1 indicated in that entry, for which indeed the sbottom 
cannot be made light but the stau can; at the smaller value jj/Mq = .5 the SU(5) boundary 
conditions must and do allow the sbottom to be light, since they contain the SO(IO) boundary 
conditions as a special case. 

Recall that in the PQ- and R-symmetric SO(IO) analysis, in which M|u(5) = 0, we were 
interested in the dependence of the allowed area on Xq and Xf/X^, as depicted in Fig. 8(a). 
The analogous contour plot of the allowed SU(5) volume is shown in Fig. 8(b). One could also 
ask for the volume allowed when the SO (10) scalar mass boundary conditions are perturbed in 
the SU(5) direction to the same extent that the SO (10) Yukawa coupling boundary conditions 
are relaxed: namely, restrict |M|u(5)/Mo| < Af/A^ — 1 (and then normalize the volume by 
dividing it by Xf/Xf - 1). The answer is very simple: up to a normalization factor of order 
unity (due to the arbitrary definition of unit area and unit volume), the contour plot of this 
restricted allowed volume is similar to that of the allowed area in pure SO(IO). Wc learn 
that, as we saw in the particular example of a light right-handed neutrino, small scalar mass 
splittings and small Yukawa coupling splittings affect the allowed area to a similar extent. 

Fig. 10 shows to what extent the PQ and R symmetries may be broken. As in Fig. 8, we 
outline the ranges of fi/Mo and M1/2/M0 which lead to a nonvanishing allowed volume in the 
space of initial scalar masses and to acceptable bottom-tau unification, for different choices 
of Xg, Xf/Xf and as{mz)- The permissible ranges extend to small 11 and M1/2 even when 
a fairly large is required, because within the SU(5)-allowed volume one can find corners 
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Table 4. Tlie characteristics of Yukawa unification witli SU(5)-type scalar mass boundary 
conditions, in analogy with Table 3. 



where m~ ~ m? ~ m| ~ ~ ^1/2 ^ ^0 so Sb is quite large. 

Cosmological Bounds The stability of the lightest supersymmetric particle (LSP, which we 
denote by x) can lead to serious cosmological bounds on the parameters of the MSSM. In this 
section we discuss these bounds in the large tan P scenario. 



As a first constraint, the LSP has to be both electrically and color neutral otherwise 
it would have been found in searches for exotic isotopes. This is typically not a problem for 
us. As we have seen, in the most interesting large tan/3 scenarios, either /i or M1/2 or both 
are considerably smaller than all the other SUSY parameters. In these regions of parameter 
space the LSP can only be a neutralino or a chargino. By a numerical study of the 4x4 
neutralino mass matrix, we find that for tan/3 3> 1 and with the LEP bounds on n and M1/2 
the LSP is always a neutralino. In the limit \fi\ + IM1/2I ^ rnz, this property can be easily 
checked by performing a perturbative diagonalization of the mass matrix. 

The second constraint arises from the LSP relic mass density p^, which must not exceed 
the critical density of the universe today pc = (1.88 x lO'^^gcm"^)/;,^. We devote the rest of 
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Fig. 10. The allowed regions of PQ and R breaking at the GUT scale, in analogy with 
Fig. 7. 



this section to this issue. We base our discussion on Ref. where the lightest neutralino 
relic density was studied but without an emphasis on the parameter regions discussed in our 
paper. Recently the LSP abundance in the large tan/? scenario was also partially discussed 
in Ref. and where our analyses overlap there is qualitative agreement. 

The contribution of x to the present VlK^ (where Q = p/pc) is determined by how fast the 
LSP annihilates when it is non-relativistic. In practice, f^^^^ roughly inversely proportional 
to the annihilation cross section a^^ at a freeze-out temperature Tp 



m 



,/20il 



In our case 

the composition of x and its annihilation properties depend crucially on p and Mi/2- Thus 
the PQ and R symmetries provide once again the right language for classifying the different 
situations. 

Let us first consider the approximately PQ-symmetric and R-asymmetric scenario mz ~ 
p <C Ml/2- In this case the LSP is predominantly a Higgsino. For > m\y the annihilation 
into W pairs proceeds with full gauge strength via t-channel chargino exchange. The rate is 
easily sufficient to ensure that Qyh"^ <^ 1. For < mw, one has to rely on annihilation 
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into standard-model fermion pairs via s-channel vector boson exchange. Now the strength of 
the X couphngs to the Z plays a role, as does coannihilation both with the second-lightest 
neutralino x' and with the lightest chargino iS^ this case via s-channel W exchange). 
Define = ± where ^ are the neutral components of the Higgsino doublets. The 
mass-eigenstate fields for the LSP x ^i-iid the second-lightest neutralino x' are given by 

X = K + o{^]h_^ X' = h- + 0(^]h. (45) 

(The plus and minus signs obviously depend on our conventions.) The isospin quantum num- 
bers of and are such that the vertices Zh^h^ and Zh_h^ vanish, while Zh^h^ has 
full gauge strength. For large M1/2, Eq- ( |45| ) implies that the XX^ vertex is suppressed by 
O {m%/ iJ,Mii2^ ~ O [rnz / relative to the xx' ^ vertex. A similar discussion applies to 
the coupling XX^W^? which is not suppressed in this limit. Furthermore, the splittings m^—m^i 
and m^ — m^+ vanish like m|/Mi/2, so x' X~^ are all present just before the LSP freezes 
out, and coannihilation is important Thus for ultraheavy gauginos the self-annihilation 



rate cr(xx ~^ ff) is negligible but coannihilations o"(xx' ~^ ff) and ^(xx^ — > z/e, . . .) are im- 
portant, since the mass splittings are very small. The result^ is that Q^h'^ <^ 1 [Q. As 
Ml/2 is lowered below ~ 400 GeV, the Zxx vertex becomes important and self-annihilation 
becomes dominant, leading once again to ^x^"^ ^ intermediate values of the gaugino 

mass, f^x^^ reaches a maximum Q^ax ~ 10%. So the PQ-symmetric case does not suffer from 
an over density of LSP 's. 

As we lower M1/2 down to ~ rriz we recover the maximally symmetric case. Now the LSP 
contains sizeable portions of both /i+ and h^, and the Zxx vertex has essentially full gauge 
strength. As a result ^x^"^ is always well below 1. 

Finally, let us study the effects of raising fi to arrive at an R-symmetric and PQ-asymmetric 
scenario mz ~ M1/2 <C /i. Now the LSP is predominantly a bino (the hypercharge gaugino): 
X ~ -B. We will collectively denote the squark and slepton masses and also by the single 
mass parameter mo. In the limit fi ~ itlq 00 the LSP is totally decoupled and ^x^"^ i^ 
extremely large: all cross sections for bino annihilation vanish at least like m^/mQ or like 
l/(/i^ tan^ /3). Clearly this poses a potential problem for this scenario. For a quantitative 
study of the possible annihilation rates, it is useful to integrate out the heavy fields and 
obtain an effective Lagrangian Cxx hino B and the SU(2) gauginos Wi, J = 1, 2, 3. We 

write Cxx = Ci + C2 + C'i + C/^ in two-component notation as 

. ^^BHi{i-&^..\b. (46) 



— + ■ 

rriA 



1 1 / - 

7 -Hu H 2" l'^''^^ + -^^^^ 

tanp rrij^ ^ 



+ h.c. (47) 



^For quantitative estimates we used the formulae of Ref. where only the effect of coannihilation with 
X' was included. We expect that accounting also for coannihilation in the charged channel @ will lower the 
final value of fi^, and therefore strengthens our conclusion that there the relic LSP abundance is sufficiently 
small. 
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= ^^BHI^p(i-^ + ---]bH^ (48) 




A = ^j:^Bf{l- — + ...]Bf (49) 
U = I^BHltp(^l-^ + ...]WH^ + h.c. (50) 

where D is the covariant derivative acting to the right, Yf and are respectively the / 

sfermion's hypercharge and mass, W = WiTj and tj are the Pauh matrices. >Ci,2,4 are obtained 
by integrating out the Higgsinos and the Hd Higgs doublet, while £3 arises from integrating out 
the sfermions. The dots represent higher derivative terms. We have approximated ~ m\. 



Notice that £4 introduces also an 0{l/n ) overlap between the LSP and the W3. The overall 
contribution of £4 to the LSP annihilation amplitudes is (9(l/yU^). In fact the amplitude at 
0{l/n'^) gets contributions also from virtual bino and wino diagrams. 

Let us start from the PQ- and R-symmetric scenario and increase fi above ~ mz- Then, for 
intermediate values of /i, the LSP annihilation cross section is determined by £2 and £4, since 
all other terms above are either suppressed by a large scalar mass {m\, m^-) or by 1/ tan/3. 
Focusing on £2 and £4, we notice that the leading C(l//x^) contribution to the amplitude is 
only in the p-wave. This contribution is determined by the first term in £2, and its p-wave 
character is easily seen by using the equations of motion in the limit of massless fermions (It 
is accurate to neglect fermion masses since we suppose x to be below the top threshold in the 
R-symmetric scenario.) This fact is of some importance since the LSP's freeze out in the non- 
relativistic regime in which p-wave cross sections suffer a further suppression ~ Tp/m^. For 

< mwi the LSP's annihilate into // pairs via Z-boson exchange, at leading order from the 
first term in £2. We find that the relic density VL^h? exceeds 1 already for /i > 200 — 300 GeV 
(depending on how close is to \mz)- The situation does not improve when > mw,z 
and all the bosonic channels WW, ZZ, Zh and hh are open and dominant. There are two 
contributions to these channels: a p-wave term cTp ~ 1//^^ from the first term in £2, and 
an s-wave term o"s ~ l//i^ from the second term in £2 and from virtual gaugino diagrams 
involving the first terms in £2 and £4. (As already stated we use £2 to explicitly display the 
suppression of s-wave processes, a fact also observed in Ref. [^. For quantitative estimates. 



we calculated cxs by returning to the partial wave amplitudes given in Ref. ■) limit 
^ mzi (Tp is dominated by final states with zero helicity (longitudinal vector bosons 
and scalars). The amplitude in this case is readily given by the annihilation into Goldstone 
and Higgs bosons from the first term in £2 via the equivalence theorem. For the purpose of 
our qualitative discussion we only kept this leading term. (We expect the complete result 
not to be drastically different in the region ~ m^; in fact we checked this explicitly for 
XX ~^ WW by using the formulae in Ref. |^.) Using these estimates for ag^p we find, again, 
that VLyJi^ > 1 for /i > 250GeV. We are thus led to the interesting conclusion that we cannot 
essentially move away from the PQ- and R-symmetric scenario towards the PQ-asymmetric 
one, if all the other superpartners and the pseudoscalar Higgs are very heavy. Nonetheless the 
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moderate ~ 200 — 300 GeV scenario is interesting, since in this case the LSP could account 
for all the dark matter and give f2 = 1. 

When fi exceeds ~ 300 GeV some other particle (namely the pseudoscalar Higgs or some 
sfermion) has to be lighter in order to avoid an overdensity of LSP's. A quick inspection of 
the above effective Lagrangian shows that, by lowering or some m^, only the annihilation 
into // can be significantly affected [the amphtude into bosons from Ci is C(l/ tan/?)]. This 
process can be mediated by (I) t-channel / exchange via £3, or (II) s-channel pseudoscalar 
Higgs exchange via Ci. In case (I) the amplitude is p-wave and the effective vertex is ~ l/m'j. 

Supposing a sfermion / is relatively light, we get the following estimate for the bino-like LSP 
relic density: 

'lOOGeV\V V fxp 



where Yf is the fermion hypercharge (so Yq = 1/3) and Vf is the dimension of the / multiplet 
(so tq = 6). Thus, leaving all the other fine tunings untouched, an acceptable fl requires an 
additional tuning to at least an order of magnitude below its natural scale. The amplitude 

for case (II) behaves like l/fim\, but is in the s-wave. In this case only the final states bb 
and fr are relevant (the lighter fermions are suppressed by the small Yukawa couplings) . We 
estimate: 



/lOOGeVV/ Y / niA Y 



\ J V350GeV/ V350GeV. 

(We are being a little sloppy in the above equations, by neglecting terms in the sfermion 
and pseudoscalar Higgs propagators, but the conclusions would not loe changed much by a 
more careful computation) Notice that even though the above Q scales like mg, the result is 
comparable to that in Eq. ([^). This is partially due to the s-wave enhancement. We see once 
again that if /x is increased from its minimally-allowed value, then some mass parameter {m\ 
in this case) must be made light to meet cosmological bounds. Typically, such a requirement 
entails some further tuning of some GUT-scale parameters. 

Which annihilation channel is more likely? Assuming SO(IO) boundary conditions on the 
scalar masses. Table 3 shows that the only particles that can be made very light (i.e. com- 
parable to niz) are the pseudoscalar Higgs and the sbottom, while for SU(5)-type boundary 
conditions Table 4 adds the stau to this list. Hence efficient fermion production through 
sfermion exchange [case (I)] requires making ~ 100 GeV or ~ 200 GeV (under SU(5)- 
type conditions). Efficient fermion production through pseudoscalar Higgs exchange [case 
(II)] need not require to be quite as light (depending on /i), but recall that lowering 
TJiA also increases the need to tune B in order to generate large tan f3. So both channel are 
roughly equally unlikely. There is, of course, the possibility that the first- or second-generation 
sfermions are light — their initial values need not be related to those of the third generation, 
and their evolution is essentially decoupled from the third generation and Higgs sectors — in 
which case they could remedy the difficulties with the R-symmetric, PQ-asymmetric scenario. 

We conclude that whenever fi is small, whether the gauginos are light or heavy, the cos- 
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mological density of the LSP is well below critical. In the large ji but light gaugino case, LSP 
annihilation is unacceptably suppressed if all the other superpartner and pseudoscalar Higgs 
masses are large. The annihilation rate can be sufficient if some of those masses are lowered, 
either through fine-tuning of m^, mi or mf , or perhaps by appealing to the yet- unspecified 
first two generations of squarks and sleptons. Needless to say, the course nature has chosen 
will be definitively revealed by future measurements of the superspectrum. 

8 Conclusions 

In this paper we have studied some of the consequences of large third-generation Yukawa 
couplings in the minimal supersymmetric extension of the standard model, subject to various 
grand unification assumptions. We have focused our attention on the two cases A!^ = A^^ = \f 
and A^ = Af ~ Af imposed at the GUT scale ~ 10^^ GeV, but most of the conclusions are 
qualitatively unchanged if these conditions are instead enforced at the Planck or string scales, 
and a few general features remain even if one only assumes that A^ ~ A5 ~ A^ at some very high 
scale. For example, the need to tune some parameters to at least one part in fifty (~ mt/mb) is 
a generic consequence of LEP bounds and the structure of the MSSM Lagrangian. This is also 
perhaps the most bothersome conclusion: as long as the top and bottom Yukawa couplings 
start out comparable at the GUT scale, there is no way to both explain the top-bottom mass 
hierarchy naturally (in the technical sense) and avoid tuning m| <ti m| . But if this one bitter 
pill is grudgingly swallowed, the remaining features of the large tan /3 scenario are intriguing. 
Let us assume for the moment that there is no theoretical bias about physics at the GUT 
scale other than the existence of a GUT. Then, if we ask that an SU(3)c x U(l)em-invariant 
vacuum as well as the experimental rate of 6 ^ 57 are to be typical rather than unlikely 
outcomes of the GUT-scale parameters, and furthermore that the prediction for mb/rrir agree 
with its experimental value as extracted using QCD sum rules, then: 

(I) the Lagrangian at the GUT scale should display approximate PQ and R symmetries; 

(II) the value of the unified Yukawa couphng should either be Xq — 0.6 — 0.7 if the Yukawas 
are exactly unified (Af = Af = A^), or be Xq ^ 0.8 if the Yukawas are significantly 
split (for example 2A^ = A^ = Ac); 

(III) if the Yukawas arc exactly unified, the soft SUSY-breaking scalar masses at the GUT 
scale should be SU(5)- but not SO(10)-type, while if the Yukawas are split then they 
can be of either type; and 

(IV) threshold corrections to the SU(3)c gauge coupling at the GUT scale must be significant 
and negative relative to the SU(2) x U(1)y couplings. 

All of these features have phenomenological, testable consequences. We expect: 

(I) light charginos and neutralinos, which may furnish tantalizing signals at LEP II and 
would definitely be seen at the LHC; and large masses for the squarks and sleptons (at 
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least of the third generation), the pseudoscalar Higgs and the charged Higgs bosons [see 
Eq. (gg)]; 

(II) either a top mass between 160 and 170 GeV if the Yukawas are exactly unified, or 
rrit > 175 GeV if they are significantly split at the GUT scale; the Yukawa splitting at 
^GUT is reflected in the superspectrum, for example in a light stop, or better yet in the 
mass combination {Xt — Xb) / {Xt + Xb) defined using Appendix A, which is very sensitive 
to both \f /Xf splitting and to the departure from SO(IO) boundary conditions on the 
soft scalar masses (see the next point); 

(III) a large value for the (almost) invariant combination /4 = /i| + /i^ + |m? — — ^mf + 
^m~, at least if the Yukawas are exactly unified, since SO(10)-type boundary conditions 
have vanishing Mlu^g-, = and 

(IV) a value of as{mz) ~ .115 somewhat below the central gauge unification prediction. 

In particular, if feature (I) is actually borne out by future discoveries, that is, if m| <^ m| 
(where ms is a typical soft-breaking scalar mass) is transformed from an unnatural assumption 
to simply an experimental fact, then the large tan/3 scenario is as natural as the small tan/3 
conventional one. The two offer very different explanations of the top-bottom mass hierarchy. 
But the large tan /? scenario offers a more robust test of the bottom-tau unification hypothesis, 
since if and are 0{1) at the GUT scale they should be much less subject to perturbations 
by other operators. In other words, the conventional scenario suffers from uncertainties in 
bottom-tau unification from physics at remote scales, while with large tan (3 the uncertainties 
are at low energies and hence are imminently accessible. As a result, only for large tan/? is 
there a tight relationship between mj, m;,, irir and the superspectrum. 

It is important to note that some of the above predictions strongly depend on the allowed 
range for the bottom mass. The uncertainty in this mass is dominated by our estimate of 
the theoretical error in the QCD sum rule extraction. If, for example, we would know that 
mb{mb) < 4.15 GeV, then the top mass would be at least 170 GeV for exact Yukawa unification 
or at least 180 GeV for split Yukawas, and the former would be disfavored because of the tuning 
mandated by its large A^. The upper bound on ^^(m^) would also be strengthened, and would 
rely less on fine-tuning arguments. Of course, the experimental uncertainty on ^^(m^) must be 
reduced to convincingly test these predictions. But it is conceivable that in the next decade 
we will know the superspectrum well enough to calculate Sb and the logarithmic threshold 
corrections; if we can also extract the bottom mass to 0{al) and measure the top mass to 
within a GeV, then precision tests of Yukawa unification at the GUT scale would be within 
our reach. 

After this work was essentially completed, the CDF and D/) collaborations at Fermilab 
announced the long-awaited discovery of the top quark |2^. The PQ- and R-symmetric sce- 
nario we have advocated predicts a top mass that agrees very well with the values determined 
by these experiments: whereas we predict 160 GeV < rrit < 190 GeV (using the approximate 
fixed-point value as an upper bound), CDF measures rrit = 176 ± 8 ± 10 GeV, while makes 
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the less precise determination rrit = 1991*121 i 22 GeV. The argument can of course be re- 
versed: the measured value of the top mass lends further support to a PQ- and R-symmetric 
Lagrangian. As the uncertainty in rrit is reduced, it will serve as an increasingly powerful test 
to distinguish the various scenarios we have considered. 

How did we arrive at the most likely set of parameters? As sketched in Table 2, if the 
hierarchy Vu ^ vd and the suppression of r(6 — > 57) are to be obtained most naturally, 
the PQ and R symmetries should be approximately valid but without making the typical 
superpartners too heavy nor the charginos and neutralinos too light: the most desirable 
superspectrum hierarchy is /i/mo ~ m^/mo ~ 1/7. (In highly-focused situations, such as 
SO(IO) with = > 1, there can be two mass hierarchies, but since such situations are 
always more fine-tuned they are not presently relevant.) From Fig. 3 and Table 1, we learn 
that the resulting value of 6f, (~ 5%) is compatible with bottom-tau unification only if: (la) 
Ag > 0.6 and mt > 160 GeV for a, ~ 0.115 or (lb) A^ > 1 and mt > 180 GeV for a, ~ 0.125, 
if the Yukawas are exactly unified; or (2) Ag ^ 0.8, rrit ^ 175 GeV and ~ 0.115 if the 
Yukawas are significantly split at the GUT scale. Then turning to Fig. 8, we conclude that 
case (la) can be saved from further tuning by allowing SU(5) scalar mass boundary conditions 
and keeping Ag below roughly 0.7, but case (lb) would always require large tuning because 
of its large Ag. The approximately-unified case (2), on the other hand, can be naturally 
obtained by either SU(5)- or SO(10)-type scalar mass boundary conditions. Scenarios (la) 
and (2) are therefore the two we have proposed as the most likely in the absence of more 
specific model-building biases. 

What GUT models would yield these preferred scenarios? Our original motivation for 
studying unified Yukawa couplings was provided by the minimal SO (10) scenario, in which 
both light Higgs doublets lie in the same 10 h multiplet. This is the case, for example, in the 
simplest implementation of the Dimopoulos-Wilczek missing VEV mechanism |^ for solving 
the doublet-triplet splitting problem. In such models, the soft SUSY-breaking parameters 
which remain after integrating out the heavy GUT sector can be rather constrained. This is 
indeed the case when SUSY breaking is communicated to the GUT sector only via gravitational 
interactions with a hidden sector. Consequently, the structure of the soft terms it tightly 
linked (see for example Ref. |^) to the GUT superpotential. It can then be shown that 



the only source of Higgs splitting for minimal missing VEV models is the D-terms, so the 
SO(10)-type scalar mass boundary conditions hold. Therefore, to allow the freedom of SU(5)- 
type boundary conditions favored in scenario (la) above while preserving the unified Yukawa 
relations Xf = = more general soft terms would be required. With such terms, as 
may be produced when there is moduli field dependence of the GUT superpotential couplings 
1^, it may be possible to induce additional, F-type splittings between the MSSM particles. If 



scenario (la) were supported experimentally, it would thus shed some light on the mechanism 
which breaks supersymmetry. On the other hand, when the Yukawa couplings are split at 
the GUT scale, as in the second favored scenario above, the tuning can always be made 
minimal by using the PQ and R symmetries. Split Yukawas would be completely expected in 
SU(5) models or as a consequence of string theory, but they could also arise in SO(IO) models 
when the light Higgs doublets originate in several SO(IO) multiplets. Note that even with 
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universal GUT scalar masses the fine-tuning can be minimal (that is, ~ 1/ tan/3) if \f is 
significantly greater than unity. 

If we do espouse some particular class of models, we may be willing to accept a scenario 
which requires tuning to better than one part in fifty. For example, the simplest SO(IO) 
scenario with D-terms as the only GUT source of Higgs doublet splittings requires roughly 
an extra order of magnitude in tuning. But theoretically it is appealing for its simplicity, and 
its tuning can be somewhat mitigated by the phenomenologically-favored light right-handed 



neutrino (which for large tan [3 does not impair bottom-tau unification [jT5|, |T9|, ^). Moreover, 
if a model were sufficiently predictive to specify the GUT-scale boundary conditions in terms 
of few unknown parameters or even none, then the conditions we have specified for proper 
electroweak symmetry breaking would either be fulfilled — in which case the model would not 
be fine-tuned but rather remarkably predictive — or not fulfilled, in which case it would be 
ruled out. We do not know of any such models at present; until a candidate is found, we can 
only offer arguments of naturalness to point us in the right direction. 

If we are willing to sacrifice some naturalness, then also the PQ or R symmetries may be 
relaxed. Note that relaxing both would lead to a large 5b and hence a light top, in contradiction 
with the recent measurements of CDF and D/); therefore the typical SUSY-breaking scalar 
mass must be significantly above niz- Also, without these symmetries other flavor-changing 
neutral current processes beyond b ^ s'y may be problematic. In any case, the requirement 
of proper electroweak breaking does not significantly favor one symmetry over the other. 
Cosmological upper bounds on the relic LSP density, however, favor a spectrum that is only 
PQ-symmetric over one that is only R-symmetric. This is of benefit to models (see for example 



Ref. |^8|) in which the /i problem is solved by generating /i radiatively from gaugino masses 
or A terms, typically leading to ~ aMi/2- On the other hand, with light gauginos and 
large /i the predominantly-bino LSP annihilates inefficiently and "overdoses" the universe. 
To reconcile the predicted LSP relic abundance with the measured age of the universe, one of 
the superpartners or the pseudoscalar Higgs must be tuned light, or else PQ symmetry must 
be partially restored by lowering /i down to 200 — 300 GeV. Of course, at the edge of the 
allowed range for these parameters, the LSP is a prime candidate for the dark matter. 

There are many aspects to the Yukawa-unified MSSM beyond the question of naturalness. 
In Appendix A we give the exact, semi- analytic solution to the complete 1-loop RG equations 
for the third generation and the Higgs sector. Understanding their behavior under various 



assumptions and boundary conditions was a prime topic in our study. Sees. [T^ and |7j 
addressed the relationship between these boundary conditions and the superspectrum, and 
presented the ranges in which /i and M1/2 must fall to allow proper electroweak symmetry 
breaking. The effects of a light (relative to Mqut) right-handed neutrino threshold were 
examined and found helpful to symmetry breaking, while it is known that they have negligible 
impact on 6 — r unification at large tan /3. The process b —>■ sj was reexamined, and the 
possibility of cancellations between various diagrams enhanced by large tan (3 were identified 
for the first time. Finally, various issues regarding proper electroweak breaking were raised 
and resolved: the two fiat directions which could destabilize the scalar potential, and the 
scales at which they could pose a danger; the constraints on the trilinear A parameters even 
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for third-generation sfermions due to the hierarchy <^ m|; and similar constraints on the 
trihnear /x couphngs which are often neglected. 

The implications of a hierarchy in the Higgs expectation values rather than in the third- 
generation Yukawa couplings are surprisingly rich. Many aspects of the MSSM are qualita- 
tively changed by this assumption, and the phenomenological consequences of these changes 
are clear and accessible to the next generation of accelerator (and perhaps dark matter) ex- 
periments. Therefore the large tan f3 scenario offers a qualitatively different alternative to the 
often-used small tan (3 "standard" supersymmetric model. We have used criteria of natural- 
ness to distinguish between the various options for achieving large tan (3. Admittedly, these 
criteria have also revealed that all large tan P models appear to require some fine-tuning of the 
GUT-scale parameters which may not be needed for small tan/9. In other respects, however, 
such as bottom-tau unification, Yukawa unification has distinct advantages over the conven- 
tional paradigm. And in the near future, most questions of naturalness will be replaced by 
solid experimental data, which will be the final arbiter of all tan (3 scenarios, large and small. 
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Appendix A: Solving the RG Equations 



The 1-loop RG equations for the parameters of the MSSM are recounted below. We use the 
notation ^ = —^'^'^j^ where fi is the mass scale, as well as m| = — 2m^ — ii\ — 2//^ and 
m\ — rriij + m\) — ii\ + 2fj,^ where mfj — fi^ + /i"^ and m% — fM^y + /x^ are the up- and down- 
type Higgs mass parameters in the scalar potential. The soft SUSY-breaking parameters run 
according to 
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(53) 



(54) 



X^^ml + m\ + ^,l + Al^X, + Al 
X^ = m}^ + ml + pil + Al^X, + Al 



Q 

X^ = m~ + rrif + fi^jj + Al = Xr + A 
^1 + 52 + ^3 where ^3 - -2 ' 



(55) 



and 5" 

-2ml p + m - - m~ -Km? 
satisfies the useful relation 



l^z^t^A~ 2m| + -|- m~ -|- m? — m~ and 5*1^2 = 
5" evolves according to = bigfS, and therefore 
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The gaugino masses are given by Mj = Mi/2{gf / Qa) where M1/2 and go are the unified GUT- 
scale gaugino mass gauge couphng, respectively. The // parameter in the superpotential runs 
according to 

^ ^ -^\l-^\l-\l + lgl + ^gi]\l^. (57) 



Finally, the evolution of the gauge couplings is given by -^gf — bigf where 61 = —33/5, 
b2 — —1 and 63 = 3 [note that we always use gi normalized as an SU(5) coupling], while the 
Yukawa couplings evolve according to 
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(58) 



We now present the solution to the RG equations for the dimensionful parameters in terms 
of integrals over the dimensionless ones, namely the gauge and Yukawa couplings. Notice that 
/X renormalizes multiplicatively, in fact by a factor of order unity, and that is does not enter 
into the RG equations of the other mass parameters. For this reason, we may just as well 
treat /i at the clcctrowcak scale as the fundamental parameter, and thus we will not need to 
refer to its GUT-scalc value or its RG evolution. 

The RG equations for the A parameters take the form 



d_ 



A = HA + M1/2GA , H 



6A? 



XI 
3A^ AXl 



(59) 



where Ga = [fgt + H + f T^dl + ^gt + f^?|> |^?^ + 3^?2') lal- The solution is given in 
terms of the "time" -ordered exponential of the integral of the matrix H, 



n = T(exp / Hdr' 



(60) 



which satisfies -^H — HH. It may easily be computed numerically or estimated analytically. 



For example, with a GUT scale of 2.5 x 10^^ GeV and aa 
by 



Af = Af = 0.6: 

n = 

.239 -.056 .005 

-.061 .278 -.048 

.029 -.290 .610 



Af = Af = 1.0: 

n = 

.115 -.042 .007 

-.048 .160 -.050 

.039 -.273 .430 



1/24, it is approximately given 



2Af = Af = 1.0 : 

n = 

I .103 -.041 .003 

-.073 .372 -.048 

V .034 -.267 .670 



(61) 
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The trilinear couplings are then related to their GUT-scale values and to the gaugino mass at 
the GUT scale via ^ 

A = HAg + Mi/2n r n-^GAdr' . (62) 

J Tn 



The coefficient of Aq is typically an order of magnitude smaller than that of M1/2 in the 
solutions for At and Ab, so we will often assume that At^t are essentially determined at the 
electroweak scale by the gaugino masses. In the maximally symmetric case, both M1/2 and Aq 
are negligible, while without the R symmetry the M1/2 contribution is large and small effects 
due to Aq do not alter any conclusion substantially. (Of course one could imagine a scenario 
with Ag much larger than the gaugino mass and tuned to a particular value, for example to 
force At and thus suppress the rate for b —>■ 57, but we will not pursue this further. Such 
a tuning is implicitly included in eb^s-y-) We include approximate numerical expressions for 
A at the end of this appendix. 

The RG equations ( ^31) for the soft-breaking masses are now easily solved by noting that 
on their right-hand sides there are only three homogeneous driving terms, the Xi, for seven 
mass parameters. By taking linear combinations of the seven masses, specifically 
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or 




= M fl 



(64) 



we may write the RG equation of fl in the form 

d_( X \ _ f H 
dri I ~ I 





(65) 



where vs = ("el'^ss'^s'^) ^"^^ ^ ^"^^ ^ expressed in terms of gf/Qc using 

Eqs. (^). The solution to these equations is straightforward [recalling also Eq. 
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Assuming once again unification at the above values of the GUT scale and ao-, and ignoring 
the Ag contribution, yields for the coefficient vector of Mlj^ 
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Finally, the RG equation for B, 
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—B^Hb-A + GbMi/s 



(67) 



(68) 



where Hb = (— 3A^, — 3A5, — A^) and Gb — {^gf + 3g'2)/fl'G) is solved by simply integrating 
over the gaugino and A contributions: 



B = B^ 



G 



J^HbHcIt'^ Ag + M,/2f [hbH f H-^Ga dr" + G^ dr' . 



(69) 



Under the same unification assumptions as before, we obtain for Xf = Xf — 0.6: B — Bq — 
(0.36, 0.33, 0.08)1g- 1-03 M1/2; for Af = Af = 1.0: B = S^- (0.41, 0.36, 0.11)1g- 1-25 M1/2; 

and for 2Af = Af = 1.0: B ^ Bg - (0.43, 0.28, 0.07)1^ - 1-08 M1/2. 

The various integrals involving the gauge and Yukawa couplings may be approximately 
evaluated analytically, since the evolution of the Qi is known and simple while the Aj may 
be approximated in various ways, in particular near the fixed-point regime. However, for 
our purposes the semi-analytic forms presented above are sufficient. To get a feel for the 
results, we can evaluate the integrals numerically. Using the same unification scale and gauge 
couplings as above, inserting the initial conditions dictated by SU(5) symmetry, and setting 
At,b,T = Ag at the GUT scale but neglecting (for ease of presentation) the small contributions 
of Ag to the soft-breaking masses, we obtain the following (approximate) explicit solutions: 





A? = Af = 


0.6 : 








m| = 


5.45 M^/2 




+ 1-41 M|63 


+ 0.38 M| 


-0.11Ms%5)-2/x2 


m\ = 


-5.23 M^/^ 


+ 1.25Ml^^ 


- 1-50 M^g^ 




+ 0.05M|u(5) + 2//2 


ml = 


4.77 Mi% 




+ 0.53Mle^ 


+ 0.09 M|. 


+ 0.01M|u(5) 


ml = 


4.87 Mi% 


-0.21 Ml, 




- 0.29 


-0.51M|u(5) 


^1 = 


5.28Mi^/2 


-0.22Ml^^ 


+ 0.56Mi'63 


+ 0.10M|. 


+ 0.23M|u(5) 




0.15Mi^/2 




+ 0.66Mi'63 


+ 0.12M|. 


+ 0.36M|u(5) 
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0.53Mi2/2 



0.08 Mlo„ + 0.83 Mle^ 



0.31 



0.64M|u(5) 



0.19 Ag + 2.2 Ml/2 
0.17 Ag + 2.3 Ml/2 
0.35 Ag- 0.13 Ml/2 
Bg - 0.76 - 1.0 Ml/2 
0.65 /x' 



G . 



Xf = 1.0 : 



5.87 
-5.61 M 



1.20 Mi^Oh + 1-61 ^163 + 0.38 M|- 



0.05M|u(5) 



1/2 

4.63 M2/2 
4.75 



1.14 M- 



lOf 



1.72 M 



163 



- 2fi' 
SU(5) + 2/^' 



+ 0.06 Ml ' ^--^ 



- 0.27 M'lo^ + 0.45 M^g^ + 0.09 M|. - 0.05 M|u(5) 



5.14 M- 



1/2 



0.17 Mi% 
0.54 M2/2 



0.23 M^o^ + 0.54 M^g^ - 0.29 M| - 0.47 M 
0.25 Mi^o^ + 0.50 MIq^ + 0.10 M^ + 0.23 M. 
0.23 M^Q^ + 0.54 M^g^ + 0.12 Mj- + 0.41 M| 



SU(5) 



SU(5) 

0.12 M^o^ + 0.77 Mi^gg - 0.31 M^ - 0.62 M|u(5) 
0.08 Ag + 2.0 Ml/2 
0.06 Ag + 2.0 Ml/2 
0.20 Ag- 0.22 Ml/2 
Bg - 0.88 Ag - I.2M1/2 
0.43 /x^ 



2A 



G 



A 



1.0 



6.19 M2/2 - 1.15 M^o^ + 1.69 Af2g^ 
-5.25 M^/2 + 1.24 M^o^ - 1.53 M^g,^ 



4.53 M: 



1/2 



0.28 M- 



0.44 M- 



163 



0.38 M|. + 0.03 M|u(5) - 2//^ 



0.07 M^ +2/.^ 



0.09 Mj. - 0.04 M 



SU(5) 
SU(5) 



5.11 M^/2 - 0.18 M^Oh + 0-64 ^1^63 - 0.29 M| - 0.54 M|u(5) 



5.27 M- 



1/2 
0.13 Mi% 

0.52^^2 



0.19 M 



SU(5) 



0.34 M 



SU(5) 



0.23 M-^o^ + 0.54 M-^g^ + 0.10 Mj- 
0.15 M^o^ + 0.70 M^g^ + 0.12 M|. 
0.07 M|o^ + 0.85 M|63 - 0.31 M^ - 0.65 M|u(5) 
0.07Ag + 1.9 Ml/2 
0.25 Ag + 2.4 Ml/2 
0.44 Ag- 0.01 Afi/2 
fiG-0.78AG- 1.1 Ml/2 
0.57 /x^. 



52 



Appendix B: The flat direction 02 



In this appendix we discuss in more detail the constraints imphed by the flat direction 02 of 
Sec. As discussed in that section, at high enough scales (A > Ahigh) the direction may 
be stabilized by nonrenormalizable operators regardless of the sign of m^. At a lower scale, 
such operators are ineffective, and a negative leads to an instability unless the linear term 



L Q 



+ m?) 



need to estimate t 



rrin 



in the scalar potential is significant at such a scale. We 
:he scale Alow down to which this linear term may be ignored, and therefore 
above which m^A) > must be enforced. 

If 77^2 (A) > for all A between Mqut and ms, there is no instability. If m2(A) < for 
some A > Ahigh, the dangerous minimum in the potential is only at field values 02 Ahigh 
(by construction of Ahigh), so to see if it is a true minimum we must run to lower scales. If 
m|(A) < when we reach A = Ahigh, then the true minimum is at 02 ~ Ahigh and leads to 
unacceptable symmetry breaking. If ^2 (Ahigh) > 0, there is no dangerous minimum at that 
scale, and we should continue running to lower scales. If m^A) gets to zero at a scale Ac- 
above the scale Alow (to be determined below) so the linear term in the potential may be 
ignored, we must minimize the full 1-loop effective potential along the fiat direction |3^. At 
1-loop order we parametrize the fiat direction by a field 02 with zero anomalous dimension, 
so (Hu) = 02^^^^^ and (L) = (j)2Zj^^^'^, where Zu,l are wavefunction renormalization coefficient 
functions satisfying the RG equations din Zu,L/dln{(j)2/ A) = 'yu,L- (Notice that {Hu) 7^ (L) 
because the D-fiatness condition which determines the VEVs of and L is corrected by 1- 
loop wavefunction renormalizations.) Then in leading-ln(02/A) approximation the full 1-loop 
effective potential (neglecting the linear term) is completely determined by the RG equation 
for calculated using Appendix A: 



m^|02p + 



87r2 



6 



3XjXt + XlXr - -glMl - QglM, 
5 



)i</>2rin 


02 




A 



m^|02|2 + Am^|02|^ln 



A 



(73) 



At the scale Ac, where m\ = 0, the above potential has the well-known Coleman- Weinberg 
minimum at (02) ~ Ac. Moreover the vacuum energy at that minimum is ~ — Am2A^, 
which is parametrically much below the usual electroweak vacuum energy —0{m%/g'^). If 
Ac and therefore (02) get too low, the linear term ml l^^l in the potential dominates, and 
the unwanted minimum disappears. To get a rough estimate of when this happens, we just 
add the linear term to Eq. ([73|) and again minimize at A = Ac. We find that the dangerous 
minimum is eliminated because of the linear term when Ac ^ i^e^^'^{m\/ Am^ii/ X^- Thus, 
m\{h) must be prevented from vanishing only above scales of order 



A- 



LOW 



AmlXb 



271 fi 

a Xf) 



(74) 



where the quantities on the right hand side are evaluated at Ac and a is a combination 
of gauge and Yukawa coupling strengths. The above is just an estimate. For instance, 2- 
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loop RG effects and finite parts in the 1-loop potential modify the numerical prefactor on 
the right-hand side by 0{1). Moreover, notice that along 02 there is a hierarchy between 
the scale of SU(2) breaking (~ {Hy)) and that of color breaking (~ {Q))- Therefore we 
expect the 1-loop corrections to ml, which we haven't included in our estimate, to be of order 
aln((Q)/Ac) ~ aln(l/a). 



Appendix C: Approximations in the effective potential 

Throughout the paper our analysis has been based on the 1-loop-RG improved tree level scalar 
potential Vo(A). In this appendix we discuss the possible relevance of a more accurate treat- 
ment which would include the full 1-loop effective potential, 1-loop threshold corrections and 
2-loop running of the soft SUSY-breaking parameters. We will show that all these refinements 
cannot change the basic conclusions of our study. 

Let us first discuss the use of the 1-loop effective potential. This becomes necessary when 
the tree-level potential is almost fiat (or even unbounded) along some direction in field space. 
Then the quantum correction stabilizes the potential at large field strength. Minimization of 
the potential then yields a vacuum expectation value for that field of order the renormalization 
scale A at which the potential becomes fiat (dimensional transmutation) p5| . (Of course we 
are assuming that at very high energy scales the tree potential is bounded from below). In 
our study, as discussed in Sec. ^Aj we need only worry about the two fiat directions 0i and 02- 
Hence the full 1-loop effective potential is only relevant for those parameter ranges when ml 
or ml become very small, namely very near the = or m^ = planes in the space of GUT 
parameters. Thus a more correct and involved calculation would only change the margins of 
the allowed parameter space but could not significantly alter any conclusions. 

The remaining improvements are given by {i) GUT-scale thresholds, {ii) SUSY-scale 
thresholds, and [iii) 2-loop running. The first type of effects are model-dependent and have es- 
sentially been encompassed by our discussions of the various boundary conditions. Let us then 
turn to the effects of the superpartner thresholds on the allowed (usually triangular) regions 
of SO(IO) parameter space — the extension to more general initial parameters will be obvious. 
The first question is where to stop the running. In the plots in Fig. 7 the running has been 
stopped at A = mz, even though in the typical scenarios the superpartners decouple closer to 
1 TeV. A more exact analysis would then add 1-loop threshold corrections to the various mass 



parameters. (In the notation of Ref . , such threshold effects would appear as terms in the 
1-loop effective potential.) These corrections are roughly proportional to amshi{ms/mz), 
where ms collectively denotes the low energy values of the MSSM mass parameters and a 
is the appropriate gauge or Yukawa coupling strength. In the absence of strong premature 
focusing (for example, Xg well below 1.2 in the SO(IO) scenario), when all the superpartner 
masses are comparable, the effect on the triangle plots is just to move the triangle boundaries 
by an amount which is roughly 0[aln{ms/mz)] times the size of the triangle itself. The size 
of the allowed region in the GUT parameter space is only slightly changed, and the same 
arguments we have made can be applied to the slightly shifted GUT parameter ranges. If 
there is strong premature focusing — that is, if Xg is just below its maximal value, the al- 
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lowed region of parameter space is very small, and some particles have masses ml well below 
the rest of the superpartner masses ttIq fj — then the threshold corrections could be significant 
when amo^^ ln(mo_j:^/mz) ~ itT'o^l- But their only relevant effect would be to slightly shift the 
maximal value of Xq, so once again no conclusions are qualitatively altered. Notice that it was 
crucial to establish that the corrections are proportional to the electroweak- scale values of the 
masses, which can be much smaller that the GUT values when any focusing is relevant. The 
effects of 2-loop running, however, are in general proportional to the value of the masses at 
higher scales: we expect corrections ~ aM^ to the low energy values of the masses (recall that 
Mo is a typical soft mass at Mqut)- When A© is large and \f — \f , the low-energy masses 
in the allowed regions of parameter space are focused to small values rn^. When the focusing 
is strong enough that ~ oMq, some effects of 2-loop running are large. For example, a 
value of Xg which leads to a small but finite allowed triangle with 1-loop running could lead to 
vanishing allowed area using 2-loop running. So when there is large focusing the more exact 
upper bound on Xq could shift somewhat — but because of the fixed-point evolution of Xt at 
large Aq, the corresponding values of rrit will not change much. We therefore expect that all 
our conclusions are robust. 
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